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Preface 


This  report  summarizes  the  work  I  did  in 
identification  of  pilot  model  parameters  using  the 
Recursive  Least  Squares  algorithm.  It  is  my  sincere  hope 
that  the  results  of  this  investigation  will  contribute  to 
the  understanding  and  future  study  of  the  realtime 
identification  of  pilot  model  parameters. 

I  would  like  to  thank  my  advisor.  Major  Daniel 
Gleason,  for  his  contributions  to  this  work.  Major 
Gleason's  guidance  and  many  hours  of  help  are  greatly 
appreciated. 
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Abstract 

air  to  ground  test  technique  was  simulated  on 
the  Plight  Dynamics  Laboratory's  LAMARS  simulation  system. 
Data  was  recorded  on  longitudinal  stick  deflection  and 
longitudinal  pipper  errors.  The  data  was  used  to  identify 
pilot  model  parameters  using  the  Recursive  Least  Squares 
(RLS)  Algorithm.  Several  different  pilot  models  and 
discretization  techniques  are  used  to  determine  which 
method  is  best  suited  to  this  task.  Meal-Smith  Theory  is 
used  to  predict  a  range  of  pilot  model  parameters  to  be 
expected  from  RLS  identification.  Pilot  model  parameters 
are  identified  using  three  aircraft  with  different  time 
delays.  The  identified  pilot  model  parameters  and  pilot 
ratings  are  compared  to  see  if  a  correlation  exist.  The 
specific  values  of  pilot  model  parameters  predicted  by 
Neal-Smith  Theory  were  not  identified.  However,  trends  in 
the  pilot  model  parameters  predicted  by  Neal-Smith  Theory, 
for  aircraft  of  increasing  time  delay,  can  be  observed  in 
the  identifications.  CT 
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REALTIME  PILOT  MODEL  PARAMETER  IDENTIFICATION 

I.  Introduction 

Background 

During  the  months  of  March  and  April  1987,  a  joint 
simulation  was  conducted  by  the  German  Flight  Test 
Organization  ( DFVLR )  and  the  Air  Force  Wright  Aeronautical 
Laboratory,  Flight  Dynamics  Laboratory,  Flight  Controls 
Division  (AFWAL/FIG) .  The  objectives  of  the  joint  program 
were  to  Introduce  the  flying  qualities  personnel  and 
simulation  personnel  to  a  flight  test  method  for 
pllot/alrcraft  analysis  (1).  The  method  was  simulated  on 
the  Large  Amplitude  Multimode  Aerospace  Research  Simulator 
System  (LAMARS).  The  objective  of  the  testing  was  to 
analyze  the  pllot/alrcraft  system  in  the  ground  attack 
mode.  The  testing  technique  employed  by  DFVLR  uses  the 
Ground  Attack  Test  Equipment  (GRATE) .  The  GRATE  system 
is  a  series  of  lights  placed  in  a  pattern  on  the  ground  to 
serve  as  targets.  In  the  ground  attack  mode  the  pilot 
must  continuously  align  the  target  sight  in  his 
head-up-display  (HUD)  with  the  target  on  the  ground.  The 
lights  are  switched  forcing  the  pilot  to  react  to  plpper 
errors  in  his  HUD.  The  light  switching  and  realignment 
action  by  the  pilot  excites  the  pllot/alrcraft  system  over 
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the  wide  range  of  frequencies  necessary  for  system 
identification.  The  terrain  board  was  modified  with 
lights  similar  to  the  ones  used  by  DFVLR  in  flight  testing 
for  the  purpose  of  conducting  a  simulation  to  determine 
the  system's  suitability  for  the  analysis  of  flying 
qualities.  Two  aircraft  were  used  in  the  simulation. 

Both  DFVLR  and  AFWAL/FIGC  provided  an  aircraft  model  and  a 
test  pilot  for  the  simulation.  Both  aircraft  models  were 
linearized  transfer  function  models.  The  German  model  was 
called  Aircraft  A  and  had  dynamics  similar  to  a  German 
Alpha-Jet.  The  model  provided  by  the  flying  qualities 
group  (AFWAL/FIGC)  was  called  Aircraft  B  and  had  dynamics 
similar  to  an  F-15.  Both  aircraft  models  and  light  bank 
software  were  Integrated  into  the  LAMARS  Simulation  System 
by  the  author.  Test  runs  were  flown  at  the  target  lights 
by  each  pilot  in  each  aircraft.  Simulation  variables  and 
pilot  ratings  we re  recorded  for  each  run.  Pilot  comments 
were  recorded  and  pilot  ratings  were  given  using  the 
Cooper-Harper  rating  scale  (2). 

ELflfcl— 

The  problem  is  to  identify  a  transfer  function  model 
of  human  pilot  dynamics.  A  simplification  of  the 
Meal-Smith  pilot  model  will  be  used.  The  parameters  to  be 
identified  are  the  pilot  gain,  lead,  lag  and  time  delay. 

In  order  to  Identify  a  transfer  function  you  must  have 
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recorded  time  histories  of  the  input  and  output  of  the 
transfer  function.  The  data  collected  in  the  joint 
AFWAL/DFVLR  simulation  is  suitable  for  this  purpose.  Of 
particular  interest  to  the  identification  of  pilot  model 
parameters  are  the  recorded  time  histories  of  control 
stick  deflection  (output)  and  pipper  errors  (input).  A 
Recursive  Least  Squares  (RLS)  algorithm  can  be  used  to 
identify  pilot  model  parameters  from  the  recorded  time 
histories  of  control  stick  deflection  and  pipper  errors. 
The  more  pilot  compensation  that  is  required  the  worse 
the  Cooper-Harper  ratings  should  become.  Once  the  pilot 
model  parameters  are  Identified  they  will  be  compared  to 
Cooper-Harper  ratings  to  see  if  a  correlation  exists.  A 
correlation  with  Cooper-Harper  ratings  would  serve  as 
validation  that  the  pilot  model  parameters  identified  by 
(RLS)  are  accurate. 

Object Ivea 

The  objectives  of  this  research  are  to:  1)  assess 
the  ability  of  the  recursive  least  squares  algorithm  (RLS) 
to  identify  pilot  model  parameters  from  operating  records 
obtained  from  the  joint  AFWAL/DFVLR  simulation.  2) 

Provide  a  test  of  the  sensitivity  of  the  established 
(RLS)  algorithm  to  noise,  strategy  changes,  biases  and 
reduced  order  models.  3)  Analyze  the  effects  of 
loop-closures  and  time  delays  and  assess  the  value  of 


continued  research  using  operating  records  from  flight 
test  experiments.  4)  To  determine  if  it  is  feasible  to 
Identify  pilot  delay  as  well  as  pilot  lead,  lag  and  gain. 
It  is  also  an  objective  of  this  research  to  determine  the 
best  method  of  discretization  for  the  task  of  identifying 
pilot  model  parameters. 

Approach 

In  order  to  become  familiar  with  the  operation  of  the 
RLS  algorithm  and  to  be  able  to  apply  it  with  confidence 
to  a  model  of  unknown  form,  it  will  first  be  applied  to 
data  generated  synthetically.  An  open-loop  simulation 
will  be  performed  to  generate  synthetic  data  from  a  pilot 
model  with  known  parameters.  Once  comfortable  with 
identification  in  the  open  loop  case  an  aircraft  model 
will  be  placed  in  series  with  pilot  and  the  pilot  will 
close  the  loop.  Identification  will  then  be  accomplished 
for  the  closed  loop  case.  Biased  error  signals  and  pilot 
remnants  are  then  added  to  the  simlation  and  the 
identification  repeated.  The  results  are  compared  to  the 
results  of  the  uncorrupted  data  to  assess  their  effects  on 
the  identification. 

Meal-Smith  Theory  can  be  used  to  determine  a  pilot 
model  for  a  given  aircraft  and  bandwidth  criterion  that 
meets  specified  standards  of  performance.  The  pilot  model 
parameters  predicted  by  Neal-Smith  Theory  is  highly 


dependent  on  the  selection  of  a  bandwidth  criteria.  Pilot 
model  parameters  will  be  calculated  for  a  number  of 
bandwidths  to  determine  a  range  of  pilot  model  parameters 
to  be  expected  from  identification  of  the  operating 
records  obtained  in  simulation. 

In  examining  the  pilot  models  it  was  found  that  there 
are  several  methods  to  approximate  the  pilot  time  delay. 
Pour  different  methods  for  the  approximation  of  the  pilot 
time  delay  are  used  and  the  results  are  compared.  Several 
different  methods  are  available  for  the  discretization  of 
the  pilot  models.  Four  methods  were  used.  They  are  the 
forward  rectangular  rule,  backward  rectangular  rule, 
Tustin's  bilinear  rule,  and  the  zero  order  hold 
approximation.  Each  method  to  approximate  the  time  delay 
is  used  giving  four  different  pilot  models  to  be 
investigated.  Each  pilot  model  is  discretized  using  all 
four  methods  for  discretization  yielding  16  separate 
representations.  These  will  be  compared  to  see  which 
yields  the  most  accurate  results. 

All  identifications  will  be  done  using  the  batch 
least  squares  and  the  recursive  least  squares  features  of 
MATRIXx.  Pilot  model  parameters  identified  by  MATRIXx 
will  be  compared  to  the  pilot  comments  and  Cooper-Harper 
ratings  obtained  during  simulation  to  see  if  a  correlation 
exists. 


II.  Technical  Background 


Data-AcquiaUlon 

To  use  the  RLS  algorithm  to  identify  a  transfer 
function  it  is  necessary  to  have  time  histories  of  the 
input  and  output  of  the  transfer  function.  This  work 
attempts  to  Identify  a  transfer  function  of  the  human 
pilot  from  time  histories  of  his  input  and  output  obtained 
from  a  realtime  pllot-in-the-loop  simulation.  The  task 
simulated  was  ground  attack.  The  ground  attack  task  was 
simulated  by  an  array  of  lights  situated  in  a  8-pin 
pattern  on  the  ground.  The  array  of  lights  is  shown  in 
Figure  1.  The  lights  are  numbered  from  one  to  eight. 

The  arrow  in  Figure  1  indicates  the  switching  of  the 
lights.  In  order  to  obtain  an  accurate  identification  of 
a  transfer  function,  the  transfer  function  must  be 
persistently  excited  over  a  wide  range  of  frequencies. 

This  is  a  well  known  fact  from  system  identification 
theory.  This  persistent  excitation  is  accomplished  by  the 
array  of  lights.  When  in  ground  attack,  the  the  test 
pilot  attempts  to  align  the  pipper  in  his  Head  Up  Display 
(HUD)  with  the  lamp  that  is  currently  illuminated.  In 
other  words,  the  pilot  attempts  to  minimize  pipper  errors. 
After  a  specified  period  of  time  the  lamp  which  is 
currently  illuminated  is  turned  off  and  another  lamp  is 
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turned  on.  This  serves  as  a  step  input  in  pipper  error  to 
the  pilot  aircraft  system.  The  lights  are  switched  at 
specified  intervals  in  order  to  obtain  the  persistent 
excitation  over  a  wide  range  of  frequencies  required  for 
identification.  For  identification  of  the  human  pilot 
model,  the  input  is  the  longitudinal  pipper  error  and  the 
output  is  the  longitudinal  stick  deflection.  Both 
variables  were  recorded  during  the  realtime  simulations 
and  will  be  used  for  this  purpose. 

Much  work  was  done  by  the  author  in  preparing  the 
simulation  for  use  at  the  Flight  Dynamics  Lab,  Flight 
Controls  Division,  Simulation  Integration  Branch,  LAMARS 
Simulation  Facility.  For  the  purpose  of  identification  of 
pilot  model  parameters  it  is  necessary  to  have  accurate 
time  histories  of  the  input  (pipper  errors)  and  output 
(stick  deflections).  However,  so  that  the  reader  will 
better  understand  how  the  data  was  collected  I  will 
present  an  overview  of  the  simulation.  Two  areas  that  are 
of  particular  interest  to  the  collection  of  data  for  the 
purpose  of  identification  of  pilot  model  parameters  are 
the  hud  and  the  aircraft  model. 

The  aircraft  model  used  in  the  simulation  was  a 
fighter  aircraft  with  dynamics  similar  to  that  of  the 
F-15.  The  transfer  functions  for  the  longitudinal 
short  period  dynamics  are  as  follows: 


9(S) 
6  (S) 


(1) 


S(S2+  4.2S  +9.0) 

oifil  =  .46  (  .0435  +  1.935)  ( 2 ) 

6( S)  (S2+  4.2S  +  9.0) 

It  is  generally  recognized  that  an  aircraft  will  receive 
progressively  lower  pilot  ratings  as  its  response  to  pilot 
input  is  delayed.  To  simulate  aircraft  of  different 
flying  qualities  a  time  delay  was  used.  The  transfer 
functions  for  the  time  delays  were  developed  by  cascading 
a  first  order  Pad*  approximation.  The  Pad*  approximation 
for  the  time  delay  is  given  as 
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The  Pad*  approximations  for  different  time  delays  are 
given  as  follows: 
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Equations  1  and  2  were  multiplied  by  Equations  4,  5  and  6. 
These  equations  together  with  the  undelayed  equations 
give  four  different  aircraft  to  be  evaluated.  During 
testing  the  pilot  was  asked  to  fly  each  aircraft  in  the 
ground  attack  task.  During  each  run,  pipper  errors  and 
stick  deflections  were  recorded.  At  the  end  of  several 
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Figure  2  Cooper-Harper  Pilot  Opinion  Rating  Scale 


runs  for  which  the  pilot  flew  the  same  aircraft  he  was 
asked  to  evaluate  it.  The  evaluations  were  made  using 
the  Cooper-Harper  rating  scale  (Figure  2).  In 
all,  71  runs  were  made  and  time  histories  and  pilot 
comments  were  recorded. 

The  HUD  used  in  the  simulation  is  representative  of 
the  HUD  used  in  the  German  Alpha-Jet  (Figure  3).  The  major 
components  of  the  display  are  airspeed,  azimuth  angle, 
radar  height,  and  pitch  angle.  Of  particular  importance 
to  the  ground  attack  task  is  the  pipper  centered  between 
the  -10  degree  and  -15  degree  pitch  markings.  The  pipper 
symbol  changes  depending  upon  which  stage  of  the  attack 
you  are  in.  The  appearance  of  the  pipper  during  each 
stage  of  the  attack  is  shown  in  Figure  4.  The  symbols  are 
displayed  as  follows  shown  in  Table  1. 


Symbol 

1 

in  distance  (X  >  1900m) 

Symbol 

2 

in  range  (X>1900m) 

Symbol 

3 

in  firing  range  (X*17S0m) 

Symbol 

4 

at  the  end  (X*750m) 

the  cross  flashes  on/off 

every  1/3  sec. 

Pipper  Symbology 


Table  1 


azimuth  angle 
(10  deg) 


airspeed 
(kt)  , 


pitch  angle 


Head  Up  Display 
Figure  3 
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The  pipper  is  fixed  in  the  lateral  direction, 
however,  it  way  move  in  the  longitudinal  direction.  Two 


effects  are  taken  into  account  when  calculating  the 
longitudinal  pipper  location.  The  first  effect  is  caused 
by  a  change  in  target  location  in  the  HUD  as  the  aircraft 
approaches  it  on  a  fixed  trajectory.  This  effect  is 
called  trajectory  shift  and  is  readily  calculated  from  the 
geometry  of  the  attack.  A  second  but  less  significant 
effect  is  the  effect  of  drag  and  gravity  on  the  bullets. 

This  effect  is  called  gravity  drop.  Even  though  these 
effects  were  taken  into  account,  the  pipper  usually  remained 
in  the  upper  half  of  the  HUD  and  did  not  move  more  than 
approximately  75  mi 11 1-radians. 

All  data  collected  during  the  simulation  effort  was 
stored  in  standard  AFWAL/FIGD  binary  format.  It  was 
required  to  convert  the  data  to  ASCII  format  and  put  in  a 
form  that  can  be  loaded  into  MATRIXx. 

The.Pllot  JlQdel 

The  pilot  model  under  study  is  a  simplification  of 
the  classical  Neal-Smlth  model.  The  basic  definitions  of 
the  pilot  model  parameters  are: 

(7) 

,  in 

,  mill-radians 


G(8) 

6 

e 


«11  .  „  ,xp"T3S  IlHi 

e(S)  *p#xp  t28+1 

longitudinal  stick  deflection 
longitudinal  pipper  error 
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t 


in 


Kp  »  pilot  gain 

Tj  »  pilot  lead 

tj  ■  pilot  lag 

t3  *  pilot  delay 


aecs 

secs 

sees 


The  objective  of  this  work  is  to  identify  the  pilot  Model 
parameters  K^,  t ^  from  the  time  histories  of  6  and 

e  obtained  from  the  simulation.  The  pilot  model  given  in 
Equation  1  is  acting  as  part  of  a  feedback  control 
system  (Figure  5).  There  are  some  basic  assumptions  involved 
in  using  the  above  pilot  model.  It  is  assumed  the  pilot 
behaves  like  a  good  servo-controller.  That  is  to  say  he 
provides  a  specified  command-response  relationship.  It  is 
assumed  that  the  pilot  suppresses  unwanted  Inputs  and 
disturbances  and  concentrates  solely  on  the  desired 
control  of  the  aircraft.  It  is  also  assumed  that  he 
reduces  the  effects  of  variations  and  uncertainties  in 
elements  of  the  control  loop  (4). 

The  describing  function  of  the  human  pilot  is  a 
linear  model.  This  a  good  assumption  for  short  periods  of 
time.  However,  the  human  controller  is  very  non-linear. 

The  non-linear  portion  of  the  human  pilot  is  called  the 
pilot  remnant  and  is  the  difference  between  the  output  of 
the  describing  function  and  actual  pilot  output.  The 
pilot  remnant  will  not  be  identified  in  this  work.  It 
will  be  assumed  however  that  any  differences  between  the 
output  of  the  pilot  model  and  the  actual  time  histories 
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Closed  Loop  System 


recorded  in  simulation  will  be  attributable,  at  least  in 
part,  to  the  unmodeled  pilot  remnant. 

There  are  several  restrictions  on  the  use  of  Equation  7 
for  pilot  modeling.  The  input  to  the  pilot  must  be 
random.  The  pilot  should  not  become  sufficiently  familiar 
with  the  task  so  as  to  be  able  to  anticipate  the  forcing 
function.  He  must  react  to  the  forcing  function  with  no 
prior  knowledge  of  what  the  function  will  do.  This  was 
accomplished  in  simulation  by  switching  the  lights  at 
different  time  Intervals  At.  Where 

2.25  <  At  <  3.15 

Different  sequences  of  lights  were  also  used.  Because  the 
sequence  of  lights  and  the  time  Interval  between  switching 
was  varied  the  pilot  was  unable  to  become  familiar  with 
the  sequence.  He  was  not  able  to  anticipate  which  light 
would  turn  on  and  when. 

The  input  to  the  pilot  is  an  error  signal.  In  the 
case  of  ground  attack  the  error  is  displayed  on  the  HUD. 
The  primary  input  to  the  pilot  is  the  longitudinal 
distance  between  the  center  of  the  aiming  reticle  and  the 
target  light  as  seen  through  the  HUD.  This  is  the  pipper 
error.  The  system  is  compensatory  in  that  only  the  error 
and  not  the  error  rate  is  displayed.  The  pilot  attempts 
to  minimize  the  error  in  minimum  time  through  control 
stick  deflections. 


Pilot  compensation  is  very  sensitive  to  aircraft 
dynamics.  A  snail  perturbation  model  which  is  linearized 
about  a  nominal  operating  point  should  be  used.  It  is 
essential  that  the  pilots  attention  be  focused  on  his 
primary  task,  minimization  of  plpper  errors  in  minimum 
tine.  His  attention  should  not  be  diverted  to  side  tasks. 

The  pilot  should  be  a  trained  and  motivated  operator.  He 
should  be  familiar  with  the  dynamics  of  the  aircraft  and 
the  task  to  be  accomplished.  He  should  be  motivated  in 
the  sense  that  his  maximum  effort  is  being  given  to  the 
task.  Should  his  attention  be  diverted  or  he  gives  less 
than  his  maximum  effort  to  the  task,  the  data  he  generates 
will  be  poor  and  the  identification  of  his  dynamics  as 
represented  by  the  pilot  model  parameters  will  be 
questionable. 

Time  Delay 

Prom  Equation  7  it  can  be  seen  that  the  pilot  model 
contains  an  exponential  term. 

“T3S  (8) 

exp 

This  term  models  the  pilot  delay  due  to  neuro-muscular 
lag.  It  is  also  called  the  transportation  lag  or  dead 
time.  In  order  to  identify  the  pilot  model  parameters  of 
Equation  7  using  the  algorithm  RLS  it  is  necessary  to 
discretize  the  equation.  This  will  be  shown  in  more 
detail  in  a  later  section.  In  discretizing  Equation  7  it 
was  found  that  the  transportation  lag  given  by  Equation  8 


was  not  easily  discretized.  It  was  necessary  to  use  an 
approximation  to  the  pilot  delay. 


8everal  approximations  to  Equation  8  will  be  used. 

The  Pad*  approximation  to  the  time  delay  qiven  in  Equation 

3  will  be  used.  Two  other  methods  will  also  be  used. 

These  methods  are  not  as  accurate  as  the  Pad*  approximation 

but  are  simpler  to  implement. 

If  t3  is  very  small,  then  the  pilot  delay  may  be 
approximated  by 

"t3S 

exp  *  1  -  t3s  (9) 

~t3S  1 

*  rr^-s  <10> 

Such  approximations  are  qood  if  t3  is  very  small  and, 

the  input  time  function  f(t)  to  the  pilot  delay  term 
is  a  smooth  and  continuous  one.  This  means  that  the 
second  -  and  higher  order  derivatives  of  f(t)  are 
small  (5). 

In  addition  to  the  two  representations  shown  above, 
another  representation  to  the  time  delay  was  used.  It  is 
the  well-known  Pad*  approximation  for  time  delay  and  is 
given  by 


-T 

exp 


>S 


1  -  <t3/2)S 
1  ♦  (t3/2  )S 


(11) 


It  is  not  known  which  of  these  representations  for  pilot 
time  delay  is  best  suited  for  the  task  of  pilot  model 
parameter  identification.  It  is  desired  to  try  all  three 
methods  as  well  as  neglecting  the  pilot  time  delay 


-19- 


altogether.  This  gives  rise  to  four  distinct 
representations  of  the  pilot  nod el  G(s). 


0.(8)  -  K 
i  P 


(TjS  ♦  1) 

(t2S  ♦  1) 


G2(S) 


Kp(l-T3S) 


(t1S  ♦  1) 
(t2S  ♦  1) 


Kn  <T1S  ♦  1) 

°3(S)  "  (1+r^S)  (t2S  ♦  1) 

(1  -  <t3/2)S)  (tiS  ♦  1) 

°4(S)  *  Kp  (1  ♦  (t3/2)S)  (t2S  ♦  1) 


(12) 

(13) 


(14) 


(15) 


Bode  plots  of  Magnitude  and  phase  shift  give  a  good 
indication  of  which  approximation  to  pilot  delay  Most 
closely  approxlMates  the  actual  delay.  This  will  be 
exaMined  More  closely  In  Chapter  4. 


PlacietAzfttlQD  Techniques 

In  order  to  apply  the  recursive  least  squares  (RLS) 
feature  of  HATRIXx  It  Is  necessary  to  discretize  the 
pilot  Models  of  Equations  12  through  15.  This  is 
necessary  because  the  RLS  algorithm  in  MATRIXx  identifies 
the  coefficients  of  the  discrete  transfer  function  not  the 
continuous.  In  order  to  discretize  the  equations  one  Must 
have  a  relationship  between  the  S  doMaln  and  the  Z  domain. 
There  are  a  number  of  such  relationships  available.  Each 
represents  a  discrete  approximation  to  a  continuous 


system.  Threa  different  methods  fox  discretization  are 
used.  They  are  shown  as  follows: 


S 


(16) 


S 


(17) 


3 


1 

T 


Z  ♦  1 

The  above  approxiaations  to  S  are  substituted  into 
the  pilot  aodela  of  Equations  12  through  15.  The 
discrete  tiae  variable  Z  is  equivalent  to  the  advance  of 
one  tiae  step. 


(18) 


Zy(KT)  ■  y(KT+T)  (19) 

Equations  16  through  18  are  different  ways  of 
approxiaating  the  process  of  differentiation  in  the  Z  - 
domain.  Equation  16  is  called  forward  integration  and 
is  equivalent  to  estiaiating  the  rate  by  looking  forward 
over  the  tiae  Interval. 


y(K.TtTJ.~  .jriKTl  .  y(KT)  il-ll 


(20) 


Equation  17  is  called  backward  integration  and  is 
equivalent  to  estiaatlng  the  rate  by  looking  backward  over 
the  tiae  Interval 

yimi.-.mT.-.Ti  .  y(ltT)  u_j.ii  m, 


The  approxlaation  given  by  Equation  18  is  the 
trazezoid  rule  or  Tustin's  bilinear  rule,  found  by 
approxiaating  the  average  rate  over  the  Interval.  An 


additional  net hod  is  also  used.  It  is  called  Hold 


Equivalence  or  Zero  Order  Hold.  The  concept  behind  this 
technique  is  that  the  output  of  the  discretized  transfer 
function  is  equivalent  to  the  output  of  its  corresponding 
continuous  transfer  function  with  the  exception  that  the 
output  is  constant  or  held  over  the  interval.  The 
discretized  transfer  function  is  formed  from  the 
continuous  transfer  function  from  the  following 
relationship  (6). 


Hho(Z) 


(1  -  Z_1)  Z 


HtS) 

S 


(22) 


The  script  z  denotes  Z-transformation,  which  is  the 
Z-domain  equivalent  to  Laplace  transformation. 

Bach  of  the  four  methods  of  discretization  described 
above  are  used  to  discretize  the  4  pilot  models  given  in 
Equations  12  through  15.  This  will  yield  16  different 
discrete  representations  of  Equation  7.  Each  of  these 
representations  will  be  used  to  identify  the  pilot  model 
parameters  and  a  determination  of  which  one  is  best  suited 
to  the  task  of  pilot  model  parameter  identification  will 
be  made. 

The  general  form  of  the  discretized  pilot  model  is 
shown  as  follows: 


0(Z) 


M0Z»  ♦  H1Z*“1  »...»  Hm 

d  zn  ♦  d.z"'1  ♦. . .♦  d„ 

ox  n 


(23) 
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The  ozder  of  the  numerator  and  denominator  are  m  and 


n,  respectively.  They  will  vary  depending  on  the 
discretization  technique  used.  In  all  cases, 

n  i  m 


and  the  coefficients 


(24) 


N 


o 


9 


and 

dQ,  dl#  ...  dn 

are  identified  by  the  RLS  algorithm.  Algebraic 
expressions  for  the  coefficients  are  obtained  when  the 
pilot  models  are  discretized.  In  all  cases  they  are 
functions  of  the  pilot  model  parameters. 


No  *  £<Kp'  Tl '  r7*  t3* 


(25) 


dn  *  £(Kp'  Tl'  T2'  t3> 

The  identification  involves  solving  these  equations  for 
Kp,  t1#  Tj,  Tj.  in  some  cases  these  equations  are 
non-linear  and  do  not  yield  unique  solutions.  These 
equations  are  developed  in  more  detail  in  Chapter  4. 


Least  Squares 

Both  Batch  Least  Squares  (BLS )  and  the  Recursive 
Least  8quares  (RLS)  are  used  for  identification.  Both 
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algorithms  are  well  established  and  have  been  in  use 

extensively  in  system  identification. 

Batch  least  squares  is  a  very  popular  approach 
that  goes  back  to  Gauss  and  Legendre  in  the 
early  19th  Century.  Recursive  identification 
algorithms  update  parameters  at  every  sample, 
as  opposed  to  batch  methods  which  operate  on 
an  entire  time  history  of  data  all  at  once. 

Recursive  algorithms  are  characterized  by 
finite  non-increasing  storage  requirements. 

They  are  typically  well-suited  for  real-time 
on-line  Identification  with  modest  processors  [81. 

Both  algorithms  form  least  squares  estimates  of  the 
coefficients  from  the  following  equation: 

6'  -  (XTX)_1XTY  (26) 

where 

9 •  •  (M  ,  N. _ N  .  d  .  d,,...d  ) 

o  i  m  o  l  n 

the  Least  Squares  estimate  of 
coefficients. 

X  «  Matrix  formed  from  input  and 
output  data 
Y  =  output  vector 

The  above  equation  is  a  simplification  presented  only 
for  background  information.  In  order  to  understand  and 
use  the  BLS  and  RLS  features  of  MATRIXx  and  be  confident 
of  the  results  it  returns,  it  was  necessary  to  understand 
the  algorithms  in  greater  detail.  Por  a  detailed 
derivation  and  discussion  of  the  batch  and  recursive  least 
squares  algorithms  refer  to  reference  (6). 


An  exaaple  best  serves  to  illustrate  the  use  of 


HATRlXx  for  identification  of  the  coefficients  of  a 
discrete  transfer  function.  Given  the  following  discrete 
transfer  function: 


H(Z)  * 


611 1 
e(  Z) 


This  can  be  rewritten  as: 


a  Z  +  a 
o 

2  -  b 


l 

1 


61H  =  ao2tal  2  aofal2  (27) 

e(Z)  Z-bx  Z'1  1-^Z"1 

The  difference  equation  can  then  be  written  as: 

6(Z)  (1  -  bx  Z-1)  =  e(Z)(ao+a1Z*1)  (28) 

6(K )  =  a  e(K)  ♦  a.e(K-l)  +  b.6(K-l) 

0  1  1 

The  discrete  time  histories  of  e  and  6  from  t*l  to  t=N 

would  yield  the  following  set: 

6(1)  -  aQe(l)  ♦  axe(0)  ♦  b16(0) 

6(2)  -  a1e(2)  ♦  aLe(l)  ♦  bj6(l)  (29) 

6(3)  *  a  e( 3)  ♦  a.e(2)  ♦  b.6(2) 
oil 


6(N)  -  a  e(N)  +  a.e(N-l)  +  b.6(N) 
o  1  l 

This  can  be  written  in  MATRIX^  format  as: 

Y  -  X©' 

where 

Y  -  16(1),6(2),6(3), ...,6(N> I 


(30) 

(31) 

(32) 
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and 


e( 1 )  e(0)  6(0) 

•(2)  e( 1)  6(1) 

e(3)  e(2)  6(2) 

X  -  .  .  (33) 

•  •  • 

•  •  • 

_  e(N)  e(M-l)  6(N-1) 

The  Batch  Least  Squares  solution  to  9'  nay  be  found  by 
typing  Into  MATRIXx: 

9*  «  XsY 

MATRIX^  will  Internally  solve  Equation  30. 

To  solve  the  problen  using  the  recursive  algorithm  is 
considerably  more  simple  when  using  MATRIXx.  One  need 
only  provide  the  algorithm  with  an  initial  guess  at  the 
coefficients  (aQ,  a1#  fcl ,  the  initial  covariance  which 

5 

defaults  to  10  ,  and  the  time  histories  of  input  end 
output.  In  this  case  one  would  input: 

Y  -  (6(1),  6(2),  6(3),  ...,  6(N) ) 

U  -  (  e(l),  e( 2 ) ,  e( 3) ,  ...,  e(N)l 
NUMO  ■  (  0  0  ) 

DENO  -(10) 

(MUM,  DBM,  PITERR1-RLS  (Y,  U,  NUMO,  DEMO,  P0) 

Where  NUMO,  DEMO  are  the  initial  guesses  at  the 
coefficients.  The  program  will  return; 


MUM  *  lao,  axJ 


DEN  *  (1,  bjJ 

The  FITERR  parameter  is  defined  as  the  sum  of  the  squared 
residuals  and  qives  an  indication  of  the  goodness  of  fit. 

The  preceding  was  done  for  the  discrete  transfer 
function  of  Equation  27,  but  may  be  extended  easily  to  a 
discrete  transfer  function  of  any  order. 

Meal -smith  Theory 

Neal-Smlth  Theory  is  a  means  of  predicting  what 
compensation  the  pilot  is  likely  to  apply  and  relating  the 
compensation  to  pilot  opinion.  The  pilot  is  represented 
by  Equation  7.  The  pilot  is  viewed  as  a  good 
servo-controller  that  adapts  himself  to  the  control 
system/aircraft  combination  by  providing  the  compensation 
required  to  achieve  the  desired  performance.  The  pilot  is 
trying  to  achieve  certain  performance  standards.  The  work 
done  by  Meal-Smith  examines  pilot  performance  for  the 
pitch-tracking  task.  The  pilot's  view  of  good  tracking  is 
that  he  be  able  to  "acquire  the  target  quickly  and 
predictably."  In  the  frequency  domain  this  is  equivalent 
to  minimizing  the  bandwidth.  The  pilot  is  also  trying  to 
achieve  two  other  other  things.  He  trys  to  minimize  the 
low  frequency  droop  and  at  the  same  time  he  trys  to 
minimize  the  resonant  peak.  These  parameters  are 
illustrated  in  Figure  6.  Standards  for  these  parameters 


-27- 


have  been  established  by  Neal-Smith.  Pot  Most 
configurations  the  minimum  bandwidth  is  3.5  rad/sec.  The 
low  frequency  droop  was  deteralned  to  be  -3db.  Both 
standards  of  performance  were  determined  somewhat 
arbitrarily.  In  this  work  the  minimum  bandwidth  will  be 
varied  to  examine  its  effect  on  optimal  pilot 
compensation.  When  using  Meal-Smith  theory  to  calculate 
optimal  pilot  compensation,  pilot  time  delay  is  assumed  to 
be  constant  and  equal  to  .3.  The  theory  is  then  used  to 
calculate  the  optimal  values  of  and  t2> 

In  order  to  obtain  closed-loop  characteristics  from 
the  open-loop  system  a  Nichols  chart  will  be  used.  The 
Nichols  chart  has  the  performance  standards  defined  on  it 
as  in  Figure  7.  This  provides  a  graphical  means  to 
determine  if  the  closed  loop  pilot/aircraft  system  meets 
the  desired  standards  of  performance.  The  Nichols  chart 
is  also  useful  for  determining  the  resonant  peak. 

Meal-Smith  showed  that  there  is  a  relationship 
between  pilot  compensation  and  pilot  ratings.  Closed-loop 
resonance  is  plotted  versus  total  pilot-compensation.  In 
general,  the  lower  the  closed-loop  resonance  and  the 
closer  to  zero  the  total  pilot  compensation  the  better  the 
pilot  rating. 

In  this  work  Neal-8mith  theory  will  be  used  to 
predict  what  values  of  Kp,  t2  are  to  be  expected  from 


pilot  model  Identification  of  the  operating  records 
obtained  from  the  simulation.  Optimal  pilot  models  will 
be  calculated  over  a  range  of  bandwidth  criterion  to 
determine  a  range  of  pilot  model  parameters  to  be 
expected.  Closed-loop  resonance  will  be  plotted  versus 
the  pilot  model  parameters  obtained  from  identification  to 
see  if  a  correlation  exists  with  the  pilot  ratings  (3). 


Standards  of  Performance  Defined 

Figure  7 
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III.  Closed-Loop  Simulation 


K 


Character  lgtlcg_of_thc..BLS,  Algorithm 

The  purpose  of  this  section  Is  to  Investigate  the 
characteristics  of  the  RLS  algorithm.  The  algorithm  will 
be  used  to  Identify  pilot  model  parameters  in  the  closed 
loop  case.  In  order  to  close  the  loop  around  the  pilot 
model  It  Is  necessary  to  simulate  the  aircraft  pitch 
response  to  control  stick  deflections  generated  by  the 
pilot  model.  A  stick  to  pitch  transfer  function  was 
chosen  that  is  representative  of  a  generic  fighter-type 
aircraft  in  the  ground  attack  mode.  A  discretized  input 
signal  found  useful  in  system  identification  is  then 
applied  to  the  pilot/alrcraft  system  and  time  histories  of 
input  and  output  are  recorded.  These  time  histories  are 
then  used  to  Identify  the  pilot  model  parameters.  The 
ability  of  the  algorithm  to  identify  pilot  model 
parameters  when  the  system  is  corrupted  by  noise  is  also 
Investigated. 

The  system  under  study  is  the  closed  loop  system 
shown  in  Figure  8.  The  pilot  model  actually  represents 
the  linear  portion  of  the  human  pilot.  The  additional 
term  6*  is  added  to  account  for  nonlinear  human  pilot 
effects.  The  term  is  the  pilot  remnant  and  is  the 
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Figure 


difference  between  the  output  of  the  pilot  Model  end  the 
actual  control  deflection. 

Two  pilot  Models  will  be  identified.  The  first  Model 
is  representative  of  a  pilot  using  excessive  lead 
coMpensation  and  is  iMpleMcnted  with: 

K  ■  2.0 

P 

Ti  -  •• 

Tj  -  .2 

The  second  pilot  Model  is  representative  of  a  pilot 
using  excessive  lag  conpensation  and  is  InpleMented  with: 


1.0 


.2 


T 


2 


CE 


.8 


The  aircraft  stick  to  pitch  transfer  function  used  is 
representative  of  a  OerMan  alpha-jet  and  is  given  as  follows 

.  V  ♦  bl -  (34) 

d(s)  s'*  ♦  SjS  ♦  SjS  ♦  aQ 

Where 


b  *  1.406  aA  »  .0985 

o  o 

bx  •  2.38  -  23.5 

Sj  ■  5.67 

The  RL8  identification  is  not  dependent  on  the  aircraft 
Model  used.  The  optiMUM  pilot  Model  and  the  closed  loop 
slMulation  both  use  the  P-15  Model.  The  discretized  pitch 
input  coMMand;  ®c  is  shown  in  Figure  9.  In  order  to 


TIME -'SIC 


Discrete  -  Error  Pitch  -  Attitude  Ccrrmand  Signal 

Figure  9 


iapleMnt  the  pilot  Model,  the  pitch  response  to  stick 
input  and  the  coMMand  input  had  to  be  discretized.  The 
coaattnd  input  was  easily  discretized  because  it  is  already 
in  a  discrete  form  as  can  be  seen  in  Figure  9.  The  pilot 
and  aircraft  Model  were  not  so  trivially  discretized.  The 
pilot  and  aircraft  Models  were  both  discretized  using  the 
forward  integration  rule. 

The  pilot  Model  May  be  written  in  discretized  forM 


G  (2) 
pl 


Where 


611 1  m  Al2  ♦  Ac 

e(Z)  Z  ♦  B, 


(r./T) 


(Tj^/T  ♦  1) 


(35) 


(K..T.  ) 


Ao  *  ’  (t2/T  ♦  1) 
(K  t./T  ♦  K  ) 


(t2/T  ♦  1) 


(35a) 


(35b) 


(35c) 


Where  AQ,  BQ  can  be  identified  by  the  RLS  algorithm. 
The  aircraft  Model  can  be  written  in  discretized  forM  as: 


OIZ) 

d(Z) 


D,Z  ♦  D 


:3  ♦  c,z2  ♦  c.z  ♦  c 

Z  1  o 


(36) 


Where 


(-3/T3  -  2a-/T2  -  a,/T) 


( 1/T3  ♦  a2/T2  ♦  a2/T  f  a0) 


( 3/T3  ♦  a,/T2 ) 


(1/T3  ♦  a2/T2  ♦  ax/T  ♦  aQ) 


(37a) 


(37b) 


V  V,  V. 


The  remnants  and  bias  were  tested  on  both  pilot 
models  providing  10  separate  identifications. 

Table  2  contains  the  results  of  the  10 
identifications. 


Data  Type 

Pirst  Pilot 

Second  Pilot  ] 

K 

T, 

T  _ 

K 

T, 

T- 

P 

1 

2 

P 

1 

2 

Actual 

2.0 

.8 

.2 

1.0 

.2 

.8 

No  Bias 

2.000 

.8001 

.2000 

.9999 

.2000 

.8000 

eBias/Const 

1.0467 

1.9199 

.2541 

1.0429 

.2264 

.9225 

eBlas/Sln 

1.3862 

1.3437 

.2346 

1.008 

.2110 

.8425 

Rem/Const 

1.5879 

1.1101 

.2209 

1.031 

.2139 

.8708 

Rem/Sin 

1.8324 

.9131 

.2097 

1.031 

.2139 

.8708 

Table  2.  Closed  Loop  Pilot  Model  Parameters  Identified 


The  second  pilot  model  did  not  seem  to  be  very 
sensitive  to  noise  and  remnants.  The  first  pilot  model 
was  more  sensitive  to  noise  and  remnants.  The 
identifications  of  r2  was  fairly  accurate  but  large 
differences  exist  in  Kp  and  t^.  In  general  the  algorithm 
seems  to  be  fairly  sensitive  to  both  the  constant  and 
sinus idal  biases  and  remnants. 


Qptl«UM_f  Uot.-Model 

There  are  two  objectives  to  this  section.  The  first 
objective  is  to  calculate  an  optimum  pilot  model  for  the 
P-15.  A  Neal -Smith  analysis  was  done  on  the  P-15  to 
obtain  the  optimum  pilot  compensation.  In  the  Neal-Smith 
analysis  a  Pad*  approximation  was  used  to  represent  the 


pilot  tine  delay.  The  optimum  pilot  model,  including  the 
Pad*  approximation  to  the  time  delay  are  discretized 
using  the  backward  integration  discretization  technique. 
The  P-15  is  discretized  by  putting  the  model  in  phase 
variable  canonical  form  and  using  the  approximation  to  the 
matrix  exponential  function.  Both  models  are  then  used  to 
simulate  the  closed  loop  pilot/aircraft  system.  The 
simulation  was  used  to  generate  time  histories  of  the 
input  and  output  of  the  optimum  pilot  model.  The  second 
objective  of  this  section  is  to  identify  pilot  model 
parameters  Including  the  time  delay.  The  batch  least 
squares  and  recursive  least  squares  algorithms  in  Matrixx 
along  with  the  time  histories  generated  in  the  simulation 
will  be  used  for  the  identification. 

The  aircraft  model  used  in  this  work  is  the  stick  to 

pitch  transfer  function  for  an  P-15  shown  below  as 

®ifi±  «  .60(1.588  ♦  1.896)  (42) 

5(S)  S(S2  ♦  4 . 28  ♦  9.0) 

The  optimal  pilot  model  is  calculated  for  the  P-15  because 
it  is  the  aircraft  used  in  the  realtime  simulation.  The 
optimal  pilot  model  parameters  may  then  be  compared  to  the 
results  obtained  from  RL8  identification  of  the  simulation 
data. 

In  order  to  calculate  an  optimum  pilot  model  you  must 


plot  the  open  loop  pilot/aircraft  transfer  function  on  a 
Nichols  chart.  In  general,  the  transfer  function  is 


where 


Kpexp 


<-.3S)  Tl8*1 

TjS  +  l 


I  _  30  ) 

and  the  Pad*  approximation  to  e  *  is  used  and  is  given 


as 


exD("*3S)  % 

*xp  *  1  ♦  .15  S 


The  first  step  is  to  plot  only  the  time  delay  multiplied 
by  the  aircraft  transfer  function  and  adjust  the  gain 
until  the  magnitude  at  «*  ■  3.0  is  approximately  -5dB,  then 
add  lead  or  lag  compensation  as  required.  In  this  section 
it  will  be  assumed  that  the  pilot  delay  is  equal  to  300 
ml 1 1 i -seconds .  The  Pad*  approximation  to  the  time  delay, 
multiplied  by  the  aircraft  transfer  function,  is  given  as 
follows 


®c(8) 


.6PU.58S  t  1,896) 
1  ♦  .158  S(S2  ♦  4.28  t  9.0) 


After  multiplying  through  by  the  Pad*  approximation  the 


transfer  function  becomes 


*<8)  B  -.1422S2*  .77748  ♦  1.376 

®c(S)  .158 4  ♦  1.63S3  ♦  5.55S2  ♦  9.08 


In  order  to  use  the  above  equation  in  a  Nichols  plot 
program  it  was  necessary  to  substitute  8>ju>  in  for  8  and 
reduce  the  above  transfer  function  to  the  following  form 


<  j«> 


A»BJ 

C+Dj 


e 


c 

where 


(.1422o>2  »  1.1376)  +  .77740)1 
<.15o>4  -  5.55«2)  +  (-1.63c*3  +  9.0«)j 

A  «  .1422c*2  ♦  1.1376 
B  «  .7774c* 

C  -  .  15c*4  -  5 . 55c*2 

D  =  -1.63c*3  ♦  9.0« 


A,  B,C,D  were  programmed  in  a  Nichols  chart  program  and  a 
plot  was  generated.  From  the  plot  it  was  evident  that 
lOdB  of  gain  was  needed  to  raise  the  plot  to  the  correct 
level  and  that  lead  compensation  would  be  needed.  A 
bandwidth  of  c*bw  *  3  was  chosen  as  the  design  point.  At 
«»BW  *  3  the  phase  angle  should  equal  -130  degrees  (3). 

But  at  <*>Bw  *  3  the  phase  is  -160  degrees.  Therefore  the 
pilot  model  should  contribute  30  degrees  of  phase  lead  at 
&>Bw  *  3.  We  have 

T^o>  »  .60  at  w  *  3.0 

*  .200  sec 

This  gives  the  following  pilot  model 

K  (t.S+1)  «  1 0 ( . 2S+1 ) 

P  * 

which  represents  pure  lead  compensation.  The  following 
represents  the  optimal  pi lot/a ircraft  system. 


gl.fi  J  . 

*c<8> 


10(.2S  ♦  1) 


1  -  .158 

.9488  «•  1.1376 

1  +  .158 

S(S2  ♦  4.28  ♦  9.0) 
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The  system  requires  further  adjustment  in  order  to  meet 
the  standards  of  performance.  After  several  Iterations  of 
adjusting  the  lead  and  gain,  the  following  system  was 
determined  to  be  optimal  and  a  Nichols  chart  is  presented 
in  Figure  10. 

«  7.4728(  .19338  ♦  1)  „±_1 , 1376 — 1 

*c  (l  ♦  . 1SSJ  [S ( S ^ ♦  4.2S  ♦  9 .0) J 

Also  shown  in  Figure  11  is  a  Bode  plot  of  the  closed  loop 

system  for  comparison  with  the  Nichols  chart.  The  plots 

of  Figure  11  show  a  system  where  the  low  frequency  drop  Is 

-3dB  with  a  closed  loop  peak  of  zero  dB.  This  is  exactly  as 

is  seen  in  the  simulation  results  that  follow.  The 

optimum  pilot  model  calculated  indicates  level  1  flying 

qualities  because  it  requires  minimum  pilot  compensation 

and  has  a  low  closed  loop  resonance. 


Plg«otlre»gn  q£  Pilot. f 


i  Time  Delay 


The  full  pilot  model  with  time  delay  is  defined  as 


Am 

e(S) 


i  -  l mis. 


T^S  ♦  1 


-pi*  <T/2)S  t2S  *•  1 
Let  T/2  equal  and  you  have 


Am 

e(S) 


1  -  t38  »  1 


P  1  *  TjS  TjS  ♦  1 

After  multiplying  out  and  rearranging  terms,  the  above 
can  be  expressed  as 


NICHOLS  CHflRT 


300j-2P0. 0-260. 0-240. 0-220. 0-200. 0-180. 0-160. 0'140. 0-120. 0-100.0  -80.0  -60.0  -40.0  -20.0  0.0 

PHRSE  ANGLE,  DEG 

Nichols  Chart  for  Optimum  Pi  lot./ Aircraft  System 


CLOSED  LOOP  COMPENSATED  SYSTEM 


(S33H030)  IdlHS  3SUHd 


o  o  o  ©  ©  O  — 

m  n  to  o  u>  ej  ®  o 


(S138I33G)  SOfUINOtM 


6  7  0  9  TOO 

FREQUENCY  ( RRD/SEC ) 


_4Ui 

•  (S) 


a„S  4  a. 8  4 


(43) 


8  4  b«  8  4  b 

l  o 


where 


-K_  r,  r. 


T2  T3 


(44a) 


K-  (t,  -  r,) 


t2  T3 


(44b) 


t2  t3 


(44c) 


T  -  4  T. 


t2  t3 


(44d) 


T2  T3 


(44e) 


equation  43  was  discretized  using  the  backward  integration 


rule  given  by  Bquation  17.  This  approximation  Cor  8  was 


substituted  into  Bquation  43.  After  multiplication  and 


rearranging  terms,  Bquation  43  becomes 


am 

e(Z) 


ApZ  4  4  A, 


(45) 


where 


Z*  4  B%Z  4  B2 


?  4  — 4  b 
mi  2  m  O 

T  T 


*2  *1 

12  ♦  .  4  ao 


(46a) 


-2a„  a. 


A1  B 


(46b) 
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(46c) 


.  fcl 

.2 


(46d) 


( 46e ) 


When  using  the  least  squares  algor itha  it  is  AQ,  A^f  A 2, 
Bj,  and  B2  that  are  identified.  The  problea  then  is  to 
relate  these  to  Kp,  r^,  t2#  t^.  Given  kQ  through  B2  froa 
the  RLS  algor i the  the  following  set  of  algebraic  equations 


relates  thea  to  the  pilot  aodel  paraaeters. 


B,  ♦  B_  -  1 


a2t2a 


A.T  ♦  2A. 


V  A  -  *2  ’ 
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I*;  -  4a  a, ) 


b.  -  b^r- 


t2  * 


In  the  calculation  o£  the  positive  square  root  was 
chosen  to  insure  that  regains  positive. 


m3  rrm + m  ^  ttutt  s  *  <  i  s  a  mbtft*  r*  ■tt#  nr*K  ^  i!  i, 


The  transfer  function  for  the  P-15  model  is  given  in 


Equation  42.  After  transforming  this  system  to  phase 


variable  canonical  form  ve  have 


1  0 


*2 1  *  0 


(47) 


0  -9  -4.2 


Bit) 


[1.1376  .948  0]  |x2 


(48) 


The  corresponding  discrete  system  is 


where 


XT(K )  *  AT  *  XT(K-l)  ♦  BT  *  UT(K-l) 

YT(K)  ■  CT  *  XT(K ) 


(49) 

(50) 


X  ♦  AT  ♦ 


(51) 
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1.  The 


MATftIXx  was  used  to  calculate  AT  and  BT  with  T  «  . 
following  Matrices  were  calculated. 

[1  .0985  .0043' 

0  -.9613  .0804 

0  -.7240  .6235 

[.00016771 
.0043  | 

.0804  J 

An  eigenvalue  analysis  was  done  on  "AT"  to  deternlne 
if  all  of  the  roots  are  inside  the  unit  circle.  The  roots 
are 

1.0 

.7924  ♦  . 1723 j 
.7924  -  .1723J 

The  root  at  1.0  was  changed  to  .9998  to  insure 
stability.  The  above  equations  were  progranMed  in  FORTRAN 
and  inplenented  in  the  s inula t ion  as  subroutine  P15. 

Plgltil.ainulatlon 

The  digital  sinulation  used  to  generate  the  tine 
histories  for  identification  is  the  sane  as  the  sinulation 
progran  used  in  the  previous  sinulation  with  the  exception 
of  the  pilot  and  aircraft  nodels.  The  progran  was 
Modified  with  the  new  pilot  nodal  that  includes  tine 
delays.  The  aircraft  nodal  was  replaced  with  subroutine 
P15.  A  listing  of  this  progran  is  contained  in  Appendix 
A.  The  sinulation  was  run  for  110  seconds  and  tabulated 


data  foe  a,  6.  9,  and  9  ia  also  containad  in 
p  c 

Appendix  A.  It  can  be  seen  froa  the  tabulated  data  that 
THBTA  approaches  THETAC  without  overshooting  as  would  be 
indicated  by  the  Bode  plot  in  Figure  11.  The  Bode  plot  of 
Figure  6  serves  to  validate  what  is  seen  in  the  data. 


r  rm  a  r  »  m  a  nrwn  *  j  *  n  re  WTrmrm 


In  order  for  the  optiaua  pilot  aodel  calculated  froa 
the  Neal-Saith  theory  to  fit  the  foraat  of  the  discretized 
general  pilot  aodel  it  was  necessary  to  add  a  lag  tera. 

The  lag  tera  that  was  added  is  (.01  •  S  ♦  1).  It  was 
deterained  that  this  tera  adds  no  aaplitude  and  negligible 
phase  lag.  It  was  included  only  so  the  fora  of  the 
optiaua  pilot  aodel  would  aatch  the  fora  of  the 
discretized  aodel.  This  would  aeke  identification  of 
pilot  lag  possible. 

The  data  generated  by  the  siaulation  is  converted  to 

a  foraet  usable  by  MATRIXx  by  prograa  CONVERT .  Given  Kp, 

Tl'  t2 *  t3'  tzom  the  optiaua  pilot  aodel,  and  the 

equations  for  the  coefficients  of  the  discretized  pilot 

aodel,  the  coefficients  aey  be  calculated  as 

Aq  ■  -.1568 
Ax  -  .5741 

Aj  -  -.3102 
Bx  -  .6909 

B2  -  -.0545 


The  MATRIXx  diary  is  contained  in  Appendix  B.  The 
coefficients  as  identified  by  MATRIX^  from  the  time 
histories  are 

A0  -  -.1569 
Ax  «  .5741 

Aj  -  -.3104 
Bl  >  .6913 

B2  -  -.0545 

This  is  a  good  identification  as  can  be  seen  by  the 
closeness  of  the  numbers.  The  fiterr  of  .2446  is  the  sum 
of  the  squared  residuals  from  the  RLS  algorithm.  A  fiterr 
of  .2446  is  very  low  by  comparison  to  much  larger  values 
of  fiterr  obtained  on  runs  where  there  was  not  a  good 
identification. 


IV.  Discretization  of  Pilot  Models 


It  is  desirable  to  determine  which  of  the 
approximations  to  the  pilot  time  delay  most  accurately 
approximates  the  actual  time  delay.  This  can  be  done  by 
comparing  the  frequency  response  of  the  pilot  model  with 
the  actual  time  delay  to  the  frequency  response  of  the 
pilot  models  with  the  approximations  to  the  time  delays. 
The  frequency  range  of  interest  for  human  pilot  dynamics 
is  from  .1  to  10  radians  per  second  (3).  In  order  to 
generate  frequency  response  plots  for  the  pilot  models  it 
is  necessary  to  have  values  for  the  pilot  model 
parameters.  The  optiaml  pilot  model  parameters  determined 
from  the  Meal-Smith  analysis  of  the  previous  section  will 
be  used  for  this  purpose.  The  optimal  pilot  model 
parameters  are 


K  -  7.47 
P 

■  .19 

Tj  ■  .01 

Tj  ■  .15 

These  values  of  the  optimal  pilot  model  parameters  are 

substituted  into  the  pilot  models  of  Bquatlon  7  and 

Equations  12  through  Equation  15.  The  resulting  pilot 

models  are  given  in  Equations  55  through  59. 

G  IS)  ■  7  47  t exo*”*38* 

O  IS)  •  7.«7  (>01g+1) 


(55) 


(56) 


7  47  ^  » 130+1 ) 
(.018+1) 


7.47  (1-.38) 


(.193+1) 

( .018+1) 


7.47 


1 

(1+.3S) 


1,118 ill 

( .018+1) 


(57) 

(58) 


n  47  (1-.15S)  ( .198+1) 

(1+.15S)  (.018+1) 


(59) 


Bode  plots  for  Gp  through  Op  were  generated  using  the 
interactive  software  package  TOTAL  (13).  A  Bode  plot  for 
Equation  55  cannot  be  generated  using  TOTAL  because  of  the 
complex  exponential  ter*  representing  the  time  delay. 
Equation  55  Must  be  solved  numerically.  In  order  to 
generate  a  solution  let  8«j«  and  substitute  into  Equation 


55. 


7.47 


1+JLSLU"! 

1+.01( j«) 


exp 


-.3  (J«> 


The  magnitude  and  phase  angle  for  the  above  expression  can 
be  given  as 


Mag  >7.47 


/ 1  t  1,13) V 
/  1  ♦  (.01)2«2 


#  ■  tan-1  (.19o)  -  tan”1  (.Olw)  -  (.3w) 

The  above  equations  for  magnitude  and  phase  angle  are  - 
programmed  in  PORTRAM  program  BASE.  The  frequency  is 
varied  from  0  to  10  radians  per  second  by  .1  radians  per 
second.  Tabular  data  for  the  magnitude  and  phase  angle 
along  with  program  BASE  are  contained  in  Appendix  B.  Bode 


plots  Cor  G_  through  G_  are  presented  in  Figure  12 
Pi  P4 

through  15  Cor  comparison. 

The  plots  oC  Figure  12  through  Figure  15  contain  the 
actual  phase  and  suignitude  calculated  from  the  above 
expressions  Cor  phase  and  Magnitude.  The  actual  phase  and 
Magnitude  can  then  be  coMpared  to  the  approxiMate  phase 
and  Magnitude  given  by  Equations  56  through  59. 

FroM  Figure  12  it  can  be  seen  that  the  Magnitude  oC 
Gp  is  identical  to  the  Magnitude  o£  Gp.  The  phase  angle 
oC  Gp  closely  approximates  the  phase  angle  of  Gp  Cor 
Crequencies  less  than  1  radian  per  second.  For 
Crequencies  greater  than  1  they  start  to  diverge.  At  a 
Crequency  oC  10  radians  per  second  there  is  approximately 
180  degrees  dlCCerence  in  the  phase  angle  of  Gp  and  Gp  . 

Gp  is  the  pilot  Model  with  the  tiMe  delay  unModeled. 
Figure  12  shows  that  Cor  unModeled  tiMe  delay  the  output 
oC  the  pilot  Model  will  lead  the  actual  output  by  180 
degrees. 

Figure  13  is  a  Bode  plot  oC  Gp  .  There  is  somb 
magnitude  distortion  at  the  higher  frequencies  as  would  be 
expected  when  adding  the  numerator  term  as  the 
approximation  to  the  delay.  At  a  frequency  of  10  radians 
per  second  the  magnitude  of  Gp  is  approximately  12db 
greater  then  the  magnitude  of  Gp.  The  phase  angle  of  Gp 
is  relatively  constant,  dropping  to  only  -15  degrees  at  10 


radians  per  second.  This  leads  the  actual  phase  angle  by 
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approximately  100  degrees  at  10  radians  per  second.  This 
is  a  large  phase  lead  but  it  is  an  improvement  over  the 


# 


180  degrees  lead  of  G_  . 

Pi 

Figure  14  is  a  Bode  plot  of  Gp  .  The  magnitude  of 

G_  is  lOdb  less  then  the  magnitude  of  at  10  radians  per 
P3  P 

second.  The  phase  angle  of  G^  is  similar  to  the  phase 

angle  of  G  and  also  leads  the  actual  phase  angle  by 
p2 

approximately  100  degrees. 

Figure  15  is  a  Bode  plot  of  Gp  .  The  magnitude  of 

G  is  exactly  the  same  as  the  magnitude  of  G_  .  This  is 
P4  P 

expected  because  the  Pads  approximation  to  the  tine  delay 

will  not  Introduce  any  distortion  in  magnitude.  The  Pade 

approximation  to  the  time  delay  will  only  introduce  a 

phase  shift.  The  phase  angle  of  Op  is  equal  to 

approximately  -55  degrees  at  a  frequency  of  10  radians  per 

second.  Op  leads  Gp  by  60  degrees  at  this  frequency. 

The  180  degrees  lead  introduced  in  0„  by  not 

pl 

modeling  the  pilot  delay,  is  the  worst  of  the  4  cases. 

G_  and  G_  provide  similar  responses  in  phase  shift  and 
P2  P3 

are  both  off  by  about  the  same  amount  in  magnitude.  The 

60  degrees  phase  lead  G_  is  the  best  of  the  4  cases. 

p4 

Since  the  responses  of  G_  and  G_  provide  equivalent 

P2  P3 

asMunts  of  accuracy,  G  is  eliminated  from  further 

p3 

consideration. 

The  step  responses  of  the  closed-loop  optimal 
pilot/aircraft  systems  are  examined  to  determine  how  well 
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the  different  discretization  techniques  approximate  the 

continuous  responses.  The  pilot  Models  of  Squat ions  56, 

57  and  59  and  the  aircraft  nodal  of  Equation  1  are  used 

to  forn  the  close  loop  system.  The  closed  loop  responses 

are  calculated  and  plotted  using  TOTAL.  Three  different 

discretization  techniques  are  used.  They  are  the  backward 

rectangular  rule,  the  Tustin's  bilinear  rule  and  the 

zero-order  hold  approximation.  G_  .  G  ,  G_  are  each 

P1  p2  p4 

discretized  using  each  method  giving  a  total  of  9  plots 
for  comparison.  A  plot  of  the  continuous  system  is  also 
presented  for  comparison  with  the  discrete  plots.  The 
plots  are  presented  in  Figures  16  through  24.  The 
sampling  rate  or  time  period  for  these  plot  is  .1  second. 

A  time  period  of  .1  was  chosen  because  this  is  the 
sampling  rate  used  in  the  realtime  simulation.  In  each 
case  the  discretization  technique  provides  a  time  response 
that  is  an  adequate  approximation  to  the  continuous.  The 
most  accurate  discretization  technique  is  the  zero-order 
hold  approximation.  The  zero-order  hold  approximation  is 
essentially  equivalent  to  the  continuous  with  the 
exception  that  the  output  is  held  constant  over  the  tine 
interval. 

In  this  section  I  will  present  the  discretization  of 


the  pilot  models  by  the  4  methods  previously  discussed. 
The  discretisation  methods  are  as  follows: 


OPTIMAL  PILOT/AIRCRAFT  SYSTEM  CONTINUOUS 


of  Optimal  Gp^s) /Aircraft  System  (continuous) 
Figure  16 


ZERO  ORDER  HOLD 


of  C^>tinal  (^.(z) /Aircraft  System 
refer  hold  approximation) 

Figure  19 
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ZERO  ORDER  HOLD 
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Forward  Rectangular  Rule 


S 


ZzX 

T 


2.  Backward  Rectangular  Rule 

ZzX 

TZ 


S 


3. 


Tustlns  Bilinear  Rule 


8 


l  ZzX 

T  Z  +  l 


4. 


Zero  Order  Hold 
(1 


Hho<Z> 


-z_1)  X 


H(8) 

8 


First  I  will  discretize  using  each  o£  the  4 


ithods . 

In  order  to  stake  the  discretization  process  easier  I  will 
put  the  pilot  smdel  in  a  convenient  form. 


G  (8) 
pi 


t  Tig  *  i 

P  T,S  ♦  1 


Vi8  »  Kn 

t28  ♦  1 


(60) 


(60a) 


Where 


.  *lg  ♦ 

bl«  ♦  bo 

(60b) 

•i  •  Vi 

(61a) 

■o-KP 

(61b) 

bl  "  t2 

(61c) 

bo  "  1 

(61d) 
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Making  the  substitution  for  S  given  by  the  forward 
rectangular  rule  into  Equation  60b  and  rearranging  terns 
leads  to  the  discrete  transfer  function  given  below  as. 


Where 


G  (Z) 
pl 


A.  Z  +  A 

_  o 

Z  ♦  B. 


(62) 


(63a) 


(63b) 

(63c) 


The  coefficients  A,,  A  .  B  are  identified  by  the  RLS 

xoo 

algorithm.  It  is  necessary  to  use  the  values  of  AQ,  AQ, 
BQ  and  Equations  61  and  63  in  order  to  solve  for  Kp, 
and  r2.  The  following  set  of  equations  solves  for  the 
pilot  model  parameters  from  Equations  61  and  63. 


*»  *  *1 
1  +  Bo 

A1T 
1  ♦  B. 


(64a) 

(64b) 


1  ♦  B. 


And 


a. 


(64c) 

(6Sa) 

(65b) 

(65c) 
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The  above  equations  are  for  the  forward  rectangular  rule. 
Making  the  substitution  for  S  given  by  the  backward 


rectangular  rule  yields, 


Kl  b,  ♦  b  T 
X  o 


-a, 


bl  ♦  boT 


-b. 


B_ 


a_  * 


bl  +  boT 

Ao  +  A1 
Bo  ♦  1 


a,  ■  - 


b,  «  - 


TAo 

V1 

TB_ 


V1 


(66a) 

(66b) 

(66c) 

(67d) 

(67e) 

(  67f ) 


The  Tustin's  bilinear  rule  yields. 

A. 


B. 


‘•l  “o* 

2b.  ♦  b  T 

1  o 

(68a) 

-  2a.  ♦  a  T 

2b.  ♦  b  T 

1  o 

(68a) 

-  2b.  ♦  b  T 

L  o 

2b,  ♦  b  T 

X  0 

(68c) 

Ag  ♦  A1 

Bo  f  1 

(69a) 

1  A1  "  Ag 

2  B_  ♦  1 

o 

(69b) 

T  <B»  -  *> 

”2  (B__  ♦  1) 

(69c) 
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Equation  sate  64,  67,  and  69  can  all  be  related  to  the 
pilot  Model  parameters  through  equation  set  65. 


For  the  zero  order  hold  approximation  the  pilot  model 
of  Equation  la  is  first  put  in  the  following  form. 


Gp  (S) 


T^S  »  1 
^P  T-S  ♦  1 


(60a) 


Vhere 


Til  [•  ♦  1/t.' 

,,(8)  -  |g  7  1/T* 


K  (8  ♦  a) 
(8  ♦  b) 


K  *=  K 

P  T2 

a  «  1/Tj^ 
b  «  1/t— 


(70a) 

(70b) 

(71a) 

(71b) 

(71c) 


The  zero  order  hold  approximation  to  Gp  is  given  by 
Equation  22. 

Vhere 


HULL  -  K(s  ♦  a) 

8  8(8  ♦  b) 

It  is  convenient  to  use  partial  fraction  expansion 

K(8+a)  lx  ♦  12- 
8(8fb)  8  S+b 


Vhere 


rl 


-  K 


A 

b 


r2 


K  (-b+a) 
*  -b 


(72) 


(73) 


(74a) 

(74b) 


The  X- teens fora  of  Equation  73  is 


H(Z) 


r1Z 

Z-l 


r^Z 


Z-r. 


(75) 


Where 


r,  -  exp 


-bT 


The  zero  order  hold  approxiaation  is  given  by 


Hho<Z>  " 


Al2  ♦ 

Z  ♦  B_ 


Where 


A1  *  rl  4  r2 
Ao  "  "r2  '  rlr3 
Bo  “  ”r3 


(74c) 

(76) 

(77a) 

(77b) 

(77c) 


The  zero  order  hold  approxiaation  was  checked  by  using 
known  values  for  Kp,  r^,  r2  and  calculating  A1#  AQ  and  BQ 
using  Equations  71,  74  and  77.  These  nuabers  were  then 
coapared  to  the  nuabers  obtained  by  using  the  discretize 
coasts  nd  found  in  MATRIX^.  The  results  are  Identical. 

This  serves  to  confira  that  the  required  Equations  71,  74 


and  77  are  correct.  Given  A  .  A.,  and  B.  froa  an 

O  L  O 

identification,  it  reaains  to  relate  these  back  to  the 
pilot  aodel  paraaeters.  The  following  set  of  equations 
will  accoaplish  this. 

Ai  ♦  \» 

-1 - *  (79a) 


1  ♦  *0 

Al»a  ’Aa 
l  ♦ 


(79b) 
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cv,v# 


(78c) 


r3  '  'Bo 


a  « 


rl  +  r2 


r1ln(r3) 

"T(rl+r2> 

ln<r3) 


T.  * 


« 


-T 

Ka/b 

1/a 

1/b 


Discretization  of  Pilot  Modal  2 


© 


G  (8) 
P2 


( r *  8  + 1 ) 
Kp(l~T38,<T2SU) 


K  t | r -8*  ♦  KJt.-t  )S  ♦  Kn 
-  D  1  3  P  1  3  p 

Tj8  ♦  1 


i,8  ♦  a.  8  ♦  a_ 

*  1  p 


»1*  ♦  bo 


Vhara 


■2 

ll 


•Kp  T1  t3 


*p  <T1  '  V 


a_  ■  K. 


b*  ■  t 


2 

1. 


77- 


(79a) 

(79b) 

(79c) 

(80a) 

(80b) 

(80c) 


(81a) 

(81b) 

(82a) 

(82b) 

(82c) 

(82d) 

(82a) 


Pilot  Model  2  was  discretized  using  the  4  Methods 
previously  discussed.  Bquation  81b  is  an  improper 
transfer  function.  After  discretization  it  is  still 
iMproper  when  using  the  forward  rectangular  rule, 
backward  rectangular  rule,  and  the  zero-order-hold 
approx lMat ion. 

When  using  Tustin's  bilinear  rule  to  discretize 
Bquation  81b  the  result  is  a  proper  discrete  transfer 
function.  Since  the  RLS  feature  of  MATRIXx  can  only  be 
used  with  proper  discrete  transfer  functions  this  will  be 
the  only  technique  used  with  pilot  Model  2.  Using  Tustins 
bilinear  rule  to  discretize  Bquation  81b  yields 


p2 


♦  A,Z  ♦  A. 


Z‘  ♦  B.Z  ♦  B 
X  o 


(83) 


Where 


A 2  * 


2b. 

T 

♦  bo 

•«. 

T  + 

*•„ 

2b. 

T  ♦ 

b. 

4a  2 

T 

2a. 

T  ♦  *o 

(84a) 


(84b) 


(84c) 
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2b 


(84d) 


’1  2b, 


-  2b, 


♦  b. 


(84e) 


The  coefficients  Aq,  B^,  Bq  are  identified  by  the 

RLS  algorithm.  The  following  set  of  equations  relates 


these  to  the  pilot  model  parameters. 


b  *  1 . 
o 


■l  -  T  [it  * 


(85a) 


(85b) 


Let  the  common  denominator  of  Equation  84a  through  84d 


equal . 


A2  ’  T  *  b0 


OS) 


Then  from  Equations  84a  through  84c  form  the  following 


itrlx  equation. 


*o  "  T  *2  Ao 


T  *2  A1 
T2a2  *2 


(87a) 


(87b) 


(87c) 


(88) 


*,v«  V*  KV  * 
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-1 


2T  4 


0  -8 


(89) 


T  -2T  4 


It  is  necesary  to  use  the  quadratic  equation  and  Equations 
82  to  solve  for  the  pilot  model  parameters. 


Tj  ■  -b  ♦ 


b  -  4ac 


The  positive  value  of  is  selected  in  each  case. 


T1  "  *  T3 


K  -  a 
P  ® 


0n  (8)  -  K 
P4  P 


f1  -  t3s1  fliULil 

|l  ♦  T38j  Lt2S  ♦  lj 


K-T.T-S*  ♦  I  (T.  -  T.)S  ♦  K. 


t2t38  ♦  (r2  ♦  Tj)8  ♦  1 


.8  ♦  a, 8  ♦  a. 


b.8  ♦  b.8  ♦  b 

*  l  o 
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(90a) 


(90b) 


(90c) 


Where 


*2  "  -Kp  T1  t3 
al  “  Kp  (T1  "  V 


•o  “  Kn 

o  p 


b2  "  T2  T3 


<t2  ♦  r3  ) 


bo  "  1 


(91c) 


(91b) 


(91c) 


(92a) 


(92b) 


(92c) 


Since  Equation  90c  is  a  proper  transfer  function  it  is 
also  a  proper  transfer  function  after  it  is  discretized. 
Equation  90c  is  discretized  using  each  of  the  4 
discretization  Methods .  After  discretization  Equation  90c 


is  of  the  following  fora. 


G_  (Z)  - 

P4 


,Z  +  A.  Z  +  A. 


Z  ♦  B.  Z  ♦  B_ 
l  o 


(93) 


The  forward  rectangular  rule  yields 


(94a) 


A1  “ 


(94b) 


-•1- 


(94c) 


♦  *1 


(94d ) 


b,  _  ^ 

.2  _  *  bo 


(94e) 


The  following  set  of  equations  solves  for  aQ,  a^  a2#  bQ, 
b^|  b2  given  ^o^  X2 • 


bo  *  1 


(95a) 


“T( 2  -  B, ) 


'1  (Bx  ♦  B2  -  1) 


(95b) 


2  (Bx  ♦  B2  -  1) 


(95c) 


Let  the  common  denominator  of  Squat ion  94a  through  94c 


equal. 


3  T2 


(96) 


•o  *  T  *3  *2 


T  A3  A1 


-  T 


-•2- 


The  following  set  of  equations  solves  for  al#  a^,  a2, 

bj/  bj  given  B^/  B^*  A^f  A2> 

b_  ■  1. 
o 

Let  a.  -  -  B  T2 

X  o 

•2  •  -  »1»2 

Then  b^,  and  t>2  can  be  found  from  the  satrix  equation. 


v  'v1'!  [M  r*i 

(B1»l)T  (B^+2)  Ibj  «j 


Let  the  cowaon  denominator  of  Equation  99a  through 


99c  equal 


The  Bquations  99a  through  99c  can  be  solved  by  the 


following  equations. 


A,  Aj  T 
A,  A,  T 

A4  \  * 


2 

2 

2 


(100) 


(101) 


(102) 


(103) 


The  Tustin's  bilinear  tula  yields 

4a,  2a, 


A1  ■ 


4b-  2b. 

— ?  +  — +  b 
T2  T  ° 


X*  2*o 

4b,  2b, 

+  *  b» 


4a,  2a. 


Ao  "  4b,  2b, 


B1  " 


2  4  — ^  4  b 
T2  T  ° 


-Ob,  2b. 


4b,  2b, 

— *  ♦  — L  ♦  b 
T2  T  ° 


4b,  2b, 


*2  "  4b,  2b, 


?  T  ^ 


(104a) 


(104b) 


(104c) 


(104d) 


(104e) 
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The  following  set  of  equations  solves  for  aQ,  a2«  bQ# 
b^,  bj  given  B^,  kQ,  A2 


bo  *  1 

•x  -  -T^Bj  ♦  2) 

e2  -  -T2(B2  ♦  1) 


Then  and  b2  can  be  found  fzoa  the  following  ssitrix 
equation. 

[2B1T  4(8^2)!  p»xl  ^  TeJ  ( 

|2(B2-1)T  4(B2*dJ  |b2J  [e2 


(105) 


bi|  .  PV  4(Br2>l  1  K] 
b2J  |_2  ( Bj-DT  4(B2+1)J  [e2J 


(106) 


Let  the  commn  denominator  of  Equation  104a  through  104c 


equal 


— ?  ♦  — *  ♦  b 
T2  T  ° 


The  equations  104a  through  104c  can  be  solved  for  aQ, 
a2#  from  the  following  set  of  equations. 


«  T4 


-  T4 


**  ** 


(107) 
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Por  Meh  of  tho  discretisation  techniques  just  presented 

V  *2'  bo#  b2'  *nd  b2  *re  aU  r*ut#d  fco  ri*  r2* 
and  tj  through  the  seat  sat  of  aquations.  Equations  91a 

through  92f  are  used  to  solve  for  the  pilot  Model 

paraaeters  once  having  solved  for  aQ,  a^,  a2,  bQ,  b^f  b 2, 

for  each  of  the  techniques.  The  quadratic  equation  is 


used  for  solving  for  r^, 


P  < 
a  ■  1. 


(109) 


then 


(110) 


■  -b  ♦ 


b  -  4ac 


the  sign  is  chosen  such  that  Is  positive. 


Tj  ■  ti  ♦  b 


t2  ■  r 


(111) 


(112) 


The  sero-order  hold  technique  is  also  used  to  discretise 
Equation  90a.  Then  using  the  sero-order  hold  technique 
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the  following  format  is  mors  convenient  then  Equation 
90c. 

I58  ♦  a. 8  ♦  a  <«| «> 

G  .  — * - 1 - SL  U13) 

P4'  1  <8  ♦ 

Where 


a_ 

IB 

Tx  T3 

(114a) 

2 

‘  “p  T2  Tj 

a. 

■ 

k  <TJ  '  lai 

(114b) 

1 

P  TjTj 

■ 

lu- 

(114c) 

*0 

T2T3 

b 

m 

1 

(114d) 

°1 

T2 

b 

m 

1 

( 114e ) 

°2 

T3 

the  xero-order  hold  approxiamtion  to  Op  is  given  by 
Equation  22  where 


Hlfll  A2g2  *  *18  *  *o  (115) 

8  "  8(8  ♦  bx)(8  ♦  b2) 


as  before  it  is  convenient  to  use  partial  fraction 


expansion. 


Where 


(US) 


(117a) 


•,bJ  -  ^b,  ♦  •„ 

r2  ■  -bj  <-bj  ♦  bj) 

(117b) 

•ib2  '  *lb2  ♦  «. 
r3  “  -bj  (-bj  4  b2) 

( 117c) 

The  X-transfora  of  equation  116  is 

r.X  e.Z  t,l 

"<*>  -  iir*  tk;*  A; 

(118) 

-b,T 

Vhoro  r4  «  oxp 

(117d) 

-b2T 

r5  -  oxp 

(117o) 

Tho  zoro-ordor  hold  approxiaation  to  0^  ia 

-  a-*-1!  *  [“s11] 

(119) 

this  can  ba  written  in  its  siaplost  fora  as 

A,*2  4  A. X  4  A 

Mho<*»  *  l - 1 - * 

X2  4  B  X  4  B 

1  o 

(120) 

Vhoro 

*2  *  *1  4  *2  *  *3 

(121a) 

tj  *  rl(-r<-rs>  ♦  ijl-l-tj)  ♦ 

(121b) 

Ao  *  *l*«rS  *  *2r»  4  ')r» 

(1210 

-  -r4-r5  (121d) 

»o  *  *4  *5  U21o) 


-if 


1 


1 


I 
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The  zero  order  hold  approx  tut  ion  for  G_  was  checked 

p4 

using  known  values  for  Kp,  T^,  T 2,  Tj  and  calculating 
Aji  A^,  Aq,  B,,  and  Bfl.  These  numbers  were  then  compared 
to  the  numbers  obtained  using  the  discretize  command 
found  in  MATRIXx.  The  results  are  identical.  This 
serves  to  confirm  that  Equations  114,  117,  and  121  are 
Indeed  correct.  Given  the  values  of  A2,  A^,  AQ,  B^,  BQ 

from  identification  using  the  RLS  algorithm  it  remains  to 
relate  these  coefficients  back  to  the  pilot  awdel 
parameters.  The  following  set  of  equations  will 
accomplish  this. 

The  quadratic  equation  may  be  used  to  solve 
equations  121d  and  121e  simultaneously  for  r^  and  r5. 

Let  a  -  1 

b  *  B, 


r5  -  -  b  ♦ 


b  -  4ac 


(122) 


r4  "  " 


>1  ’  r5 


(123) 


The  above  values  of  r^  and  r$  can  be  substituted 
into  Equations  121a  through  121c.  Equations  121a  through 


121c  can  be  written  in  matrix  form  as 


(123a) 


*  « 

11  1 

•  ■ 

rl 

k2 

-r4-r5  -l-r5  -l-r4 

r2 

m 

A1 

r4*r5  r5  r4 

»  « 

r3 

solving  for  r^,  r2  and  r3  and  three. 


5 


and  solving  equations  117  yield. 

.  ?n(r4) 


L«t 


then 


el 

•2 


-  ll><g5) 


•o  -  bl  b2  rl 


“  •©  "  b2<_b2  +  bl>r3 
"•o'  bl("bl  +  b2,r2 


now  use  equations  114  to  relate  a2,  a^,  sQ,  b^,  b2 
to  the  pilot  Model  parameters. 


(123b) 
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T1  * 


t2  " 


t3  * 


KP  ’  b,  b2 


There  are  a  total  of  nine  discretizations  done  in  this 


section.  Bach  discretization  provides  equations  for  going 


froa  the  continuous  representation  of  a  transfer  function 


to  its  discrete  representation  and  back.  For  Gp  all  4 


Methods  are  used.  Por  Op  only  Tustin's  bilinear  rule  is 


used.  Por  G  all  4  Methods  are  used,  the  coefficients 


to  the  discrete  transfer  functions  will  be  deterained  froa 


the  tlae  histories  obtained  froa  the  realtlae 


pilot-in-the-loop  simulation  discussed  in  the  technical 


background  section.  These  nuaerlcal  values  for  the 


discrete  pilot  aodel  coefficients  will  be  used  in 


conjunction  with  the  equations  developed  in  this  section 


to  deteralne  the  continuous  pilot  aodel  paraaeters  Kp, 

Tj,  and  t3.  Three  FORTRAN  programs  were  written  for  this 
purpose.  The  accuracy  of  these  prograas  were  checked  by 


calculating  the  values  of  the  discrete  coefficients  froa 


known  values  of  Kp,  t 2,  and  t3.  These  values  of  the 


discrete  coefficients  were  then  used  to  calculate  the 
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values  of  Kp,  Tj,  t2  and  r3  froa  tha  aquations  developed 
In  this  section.  In  each  case  the  equations  worked 
correctly.  In  the  cases  where  the  quadratic  equation  was 
used  to  evaluate  a  parameter  only  1  solution  to  the 
quadratic  equation  produced  the  correct  result.  The 
solution  to  the  quadratic  equation  that  produced  the 
correct  result  for  the  test  case  is  the  one  that  will  be 
used  when  the  programs  are  used  to  solve  for  pilot  model 
parameters  from  discrete  coefficients  obtained  from 
applying  the  RLS  algorithm  to  the  realtime  simulation 
data. 


V.  Neal  Smith  Theory 


Neal-Smith  Theory  is  used  to  predict  the  pilot 
compensation  required  for  a  specific  aircraft  and 
bandwidth  criterion.  The  bandwidth  criterion  selected  is 
somewhat  arbitrary.  Neal-Smith  Theory  will  be  applied  to 
the  aircraft  model  of  Equation  1  over  a  range  of 
bandwidth.  The  bandwidth  used  in  the  analysis  will  range 
from  2.5  to  4.0  radians  per  second.  This  will  provide  a 
range  of  pilot  model  parameters  to  be  expected  from  the 
identifications  obtained  from  the  realtime  simulation 
data. 

Neal-Smith  theory  is  applied  in  the  same  manner  as 
was  demonstrated  in  the  closed-loop  simulation  section. 

In  this  section  it  is  repeated  for  several  values  of 
bandwidth. 

The  pilot  model  used  is  given  in  Equation  7.  The 
aircraft  model  used  is  given  in  Equation  1.  Equation  1  is 
the  P-15  pitch  transfer  function  with  no  time  delay. 

Figure  28  is  a  Nichols  Chart  for  the  pilot/aircraft  system 
optimized  for  a  bandwidth  of  2.5  radians  per  second.  The 
circles  on  the  plot  are  spaced  at  intervals  of  .5  radians 
per  second.  The  circle  that  represents  2.5  radians  per 
second  is  as  close  as  it  will  come  to  both  standards  of 
performance.  Figure  29  is  a  Nichols  Chart  for  the 


NICHOLS  CHRRT 
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Figure  28 
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pilot/aircraft  ayatea  opt i aired  for  a  bandwidth  of  3.0 
radiana  par  aacond.  The  Meal-Salth  analyaia  la  dona  again 
for  3.5  and  4.0  radiana  par  aacond.  Tha  reeults  of  all 
four  optiaizationa  ara  praaantad  in  tabla  3. 


"hw 

*p 

T1 

t2 

2.5 

7.0323 

.1200 

.01 

3.0 

7.4728 

.1933 

.01 

3.5 

6.6464 

.3571 

.01 

4.0 

4.8417 

.6252 

.01 

Optiaal  Pilot  Modala  Calculated  Over 
a  Range  of  Bandwidth  Criterion 
Tabla  3 


VI.  Identifications 


Analysis 

The  aircraft  Model  longitudinal  dynaaics  are 
represented  by  Bquatlons  1  and  2.  The  aircrafts  tiMe 
delay  is  represented  by  Equations  4,  5  and  6.  Aircraft  of 
different  tine  delays  are  Modeled  by  Multiplying  Equations 
4,  5  or  6  by  Equations  1  and  2.  These  different  aircraft 
are  then  flown  at  the  light  bank  and  tine  histories  are 
generated  for  identification.  In  this  section,  4  different 
runs  are  exanlned  in  great  detail.  Table  4  describes  the 
aircraft  and  the  light  switching  sequence  used  for  each 
run. 


Run 

XlER  Delay  (Sec) 

Light  Sequence 
BeoMant 

12  2 

1 

.0 

4  6  8 

2 

.0 

8  6  8 

4 

.1 

4  6  8 

5 

.2 

8  6  8 

Description  of  Test  Runs 


Table  4 

Run  4  and  run  5  are  Included  for  conparison  with  run  1  and 
run  2.  The  light  sequence  for  a  particular  run  can  be  seen 
by  conparlng  the  sequences  in  Table  4  to  the  light  bank  in 
Figure  30. 
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These  stquancet  war*  chosen  because  they  approximate 
a  purely  longitudinal  task.  Data  Cor  other  sequences  Is 
available  but  contains  segments  of  purely  lateral 
switching  which  Is  not  needed  Cor  ldentl Cleat Ion  oC  the 
longitudinal  pilot  model.  The  runs  described  in  Table  4 
were  chosen  because  alrcraCt  oC  diCCerent  tine  delays  were 
Clown  at  the  same  light  sequences.  This  will  allow  direct 
comparison  oC  pilot  model  parameters  Cor  alrcraCt  oC 
diCCerent  time  delay  that  ara  perCorming  the  same  task. 
Meal-Smith  Theory  predicts  an  increase  in  pilot  lead  and  a 
decraas*  in  pilot  lag  when  the  alrcraCt  time  delay  is 
Increased. 

Runs  1  and  2  both  employ  the  alrcraCt  with  zero  time 
delay.  Runs  1  and  2  will  be  used  as  a  basis  oC  comparison 
Cor  Runs  4  and  5.  Run  4  employs  the  alrcraCt  with  100 
milli-seconds  oC  time  delay  and  Run  5  employs  the  alrcraCt 
with  200  milli-seconds  oC  time  delay.  Run  1  and  run  4 
were  conducted  using  the  same  light  sequence.  Run  2  and 
run  5  were  also  conducted  using  the  same  light  sequence. 

In  order  to  apply  the  RL8  algorithm  it  was  necessary  to 
break  the  data  up  into  3  segments  Cor  each  run.  Bach 
segment  corresponds  to  the  time  period  where  the  pilot  is 
attempting  to  minimise  the  plpper  errors  to  a  particular 
light.  Bach  time  the  light  is  switched  a  separate 
identi Cleat Ion  was  done  on  that  segment  oC  the  data.  The 


••9Mnts  ace  tabled  1,  2  and  3  as  shown  in  Table  4.  The 
Identification  was  done  in  this  Manner  so  that  the  RLS 


algor l the  would  not  be  used  on  data  with  large 
discontinuities  in  the  plpper  errors.  Breaking  up  the 
data  into  separate  segments  will  also  allow  for  direct 
coMparison  of  the  pilot  Model  parameter  identified  for  a 
particular  segaent. 

An  alternate  Method  for  taking  into  account  the  pilot 
tine  delay  will  also  be  used.  Since  pilot  Model  1 
provides  no  Means  of  identifying  the  pilot  delay,  the  data 
streaas  will  be  shifted  in  tiMe  before  application  of  the 
RLS  algor ithM  in  order  to  account  for  the  delay  in  pilot 
output.  Two  different  pilot  tiMS  delays  will  be  exaained, 
T  ■  .Is  and  T  ■  .2s.  A  pilot  delay  of  T  *  .38  was  also 
tried  but  the  fiterr  be  cams  Much  greater  and  the  values  of 
and  t2  becasw  unrealistic.  Figures  31  through  38  are 
tiMe  histories  of  the  control  stick  input  and  pipper 
errors  for  each  run.  The  longitudinal  control  stick 
deflections  are  labled  BLV  and  are  given  in  Mllli-saters. 
The  longitudinal  pipper  errors  are  labled  TPIPSRRY  and  are 
given  in  Mill 1 -radians. 

IfliPtUlcatlopa 

The  pilot  Model  paraMeters  identified  using  pilot 
Model  1  with  no  representation  for  the  t Imo  delay  are 
contained  in  table  5  through  10.  Also  contained  in  table 


TIME  SEC 

Longitudinal  Control  Stick  Deflection 


Run  2 
Figure  33 
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Run  4 


(! 


Figure  36 


S  through  10  «r«  the  pilot  aodel  parameters  identified 
when  the  data  streams  are  shifted  to  account  for  the  pilot 
delay.  Pilot  time  delays  of  .1  seconds  and  .2  seconds  are 
represented.  The  identification  using  Pilot  2  are 
contained  in  Table  11  and  the  identifications  using  Pilot 
4  are  contained  in  Table  12  and  Table  13. 

In  Table  5  through  Table  13,  the  following 
abbreviations  are  used: 

FRR  -  Forward  rectangular  rule 
BRR  -  Backward  rectangular  rule 
TUST  -  Tustin's  bilinear  rule 
ZOH  -  Zero  order  hold 


Pilot  •!  Sequent  11 


Run  91 

Run  84 

K 

II 

T. 

K 

II 

12 

Mo  data  shift 

PRR 

57.96 

.1607 

1.6778 

47.24 

.2607 

1.5924 

BRR 

57.96 

.1357 

1.6528 

47.24 

.2357 

1.5674 

TUST 

57.96 

.1482 

1.6653 

47.24 

.2482 

1.5799 

ZOH 

57.96 

.1595 

1.6653 

47.24 

.2587 

1.5798 

Data 

8hifted 

for  pilot  delay  of 

.Is 

PRR 

41.49 

.5742 

1.9084 

1.88 

6.1888 

-.6394 

BRR 

41.49 

.5742 

1.8834 

1.88 

6.1638 

-.6644 

TU8T 

41.49 

.5617 

1.8959 

1.88 

6.1763 

-.6519 

ZOH 

41.49 

.5704 

1.8959 

1.88 

6.3090 

-.6518 

Data 

shifted 

for  pilot  delay  of 

.2s 

! 

i 

i 

PRR 

22.95 

1.2017 

3.0120 

BRR 

22.95 

1.1767 

2.9870 

Mum 

rical  difficulty  i 

due 

to  the  closeness 

TUST 

22.95 

1.1892 

2.9995 

of  the  discrete  i 

coefficients. 

ZOH 

22.95 

1.1967 

2.9995 

Pilot  aodol  paraaeters  identified  for  Run  1  end  Run  4  (Segment  1) 


Table  5. 
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rtV»',V-  v 


n.i't  ftn,>  trti 


i»,V.  >'■  t 


Pilot  #1  8egnent  12 


Run  >1 


Mo  data  shift 

K 

II 

PRR 

10.11 

-.5869 

BRR 

10.11 

-.6119 

TU8T 

10.11 

-.5994 

ZOH 

10.11 

-.5673 

Data 

8hifted 

for  pilot 

delay  of 

PRR 

9.56 

-.0336 

BRR 

9.56 

-.0586 

TU8T 

9.56 

-.0461 

ZOH 

9.56 

-.0331 

Data 

shifted 

for  pilot 

delay  of 

PRR 

9.26 

.0457 

BRR 

9.26 

.0207 

TU8T 

9.26 

.0332 

ZOH 

9.26 

.0452 

ll 

K 

Run  84 

T. 

rJL 

.3771 

8.92 

-.6225 

.3264 

.3521 

8.92 

-.6475 

.3014 

.3646 

8.92 

-.6350 

.3139 

.3644 

8.92 

-.5983 

.3137 

.Is 

.9225 

11.49 

.0006 

.6281 

.8975 

11.49 

-.0244 

.6031 

.9100 

11.49 

-.0119 

.6156 

.9100 

11.49 

.0005 

.6156 

.2s 

1.1261 

10.71 

-.0100 

.4836 

1.1011 

10.71 

-.0350 

.4586 

1.1136 

10.71 

-.0225 

.4711 

1.1136 

10.71 

-.0097 

.4709 

Pilot  aodtl  par«Mt«rs  identified  fox  Run  1  and  Run  4  (Segaent  2) 

Table  6. 
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Pilot  II  SeqMnt  93 

Run  #1  Run  14 


No  data  shift 

II 

IX 

K 

lx. 

lx 

PRR 

10.8 

-.4988 

.7310 

9.40 

-.4900 

.4744 

BRR 

10.8 

-.5238 

.7060 

9.40 

-.5150 

.4494 

TUST 

10.8 

-.5113 

.7185 

9.40 

-.5025 

.4619 

XOH 

10.8 

-.4902 

.7184 

9.40 

-.4770 

.4618 

Data 

8hifted 

for  pilot 

delay  of 

.Is 

PRR 

14.57 

-.2298 

.7102 

12.97 

-.1493 

.5365 

BRR 

14.57 

-.2548 

.6852 

12.97 

-.1743 

.5115 

TUST 

14.57 

-.2423 

.6977 

12.97 

-.1618 

.5240 

ZOH 

14.57 

-.2258 

.6977 

12.97 

-.2458 

.5239 

Data 

shifted 

for  pilot 

delay  of 

.2s 

PRR 

14.79 

-.1320 

.7246 

12.53 

-.0780 

.4873 

BRR 

14.79 

-.1570 

.6996 

12.53 

-.1030 

.4623 

TUST 

14.79 

-.1445 

.7121 

12.53 

-.0905 

.4748 

ZOH 

14.79 

-.1297 

.7121 

12.53 

-.0759 

.4747 

Pilot  model  parameters  Identified  for  Run  1  and  Run  4  (Segment  3) 

Table  7. 
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Pilot  II  8ega»nt  II 


Run  12  Run  15 


Mo  data  shift 

K 

IX 

lx 

K 

IX 

lx 

FRR 

8.12 

.8209 

-.9294 

2.62 

5.6927 

-.8306 

BRR 

8.12 

.7959 

-.9544 

2.62 

5.6677 

-.8556 

TU8T 

8.12 

.8084 

-.9419 

2.62 

5.6802 

-.8431 

ZOH 

8.12 

.8319 

-.9418 

2.62 

5.7779 

-.8430 

Data 

Shifted  for 

pilot 

delay  of 

.Is 

PRR 

6.46 

1.4589 

-1.0776 

10.75 

-8.4349 

-.9294 

BRR 

6.46 

1.4339 

-1.1026 

10.75 

-8.4599 

-.9544 

TU8T 

6.46 

1.4464 

-1.0901 

10.75 

-8.4474 

-.9419 

ZOH 

6.46 

1.4758 

-1.0900 

10.75 

-8.5479 

-.9418 

Data 

shifted  for 

pilot 

delay  of 

.2s 

PRR 

7.42 

1.5991 

-1.2825 

7.94 

-3.9000 

-1.5625 

BRR 

7.42 

1.5741 

-1.3071 

7.94 

-3.9250 

-1.5875 

TU8T 

7.42 

1.5866 

-1.2946 

7.94 

-3.9125 

-1.5757 

ZOH 

7.42 

1.6147 

-1.2945 

7.94 

-3.9311 

-1.5750 

Pilot  Modol  paraMtan  identified  for  Run  2  and  Run  5  (8egaent  1) 


Table  8 


$ 

Pilot  II  S«qMRt  12 


Run  12 

Run  85 

K 

lx 

lx 

K 

ll 

lx 

Mo  data  shift 

PRR 

12.04 

.2158 

.9259 

10.75 

.0666 

.9615 

BRR 

12.04 

.1908 

.9009 

10.75 

.0416 

.8365 

TUST 

12.04 

.2033 

.9134 

10.75 

.0541 

.9490 

ZOH 

12.04 

.2129 

.9134 

10.75 

.0657 

.9490 

Data  Shifted 

for  pilot 

delay  of 

.Is 

PRR 

11.95 

.3404 

.9804 

10.65 

.2419 

1.2626 

BRR 

11.95 

.3154 

.9554 

10.65 

.2169 

1.2376 

TUST 

11.95 

.3279 

.9679 

10.65 

.2294 

1.2501 

ZOH 

11.95 

.3361 

.9678 

10.65 

.2395 

1.2501 

Data  shifted 

for  pilot 

delay  of 

.2s 

PRR 

11.44 

.4552 

1.0730 

11.15 

1.3132 

3.0488 

BRR 

11.44 

.4302 

1.0480 

11.15 

1.2882 

3.0238 

TUST 

11.44 

.4427 

1.0605 

11.15 

1.3007 

3.0363 

ZOH 

11.44 

.4499 

1.0604 

11.15 

1.3078 

3.0363 

Pilot  model  parameters  Identified  for  Run  2  and  Run  5  (Segment  2) 


Table  9 


Pilot  »1  Stgaant  13 


Mo  data  shift 

PRR 

Run  12 

lift  II 

37.4  .1382 

T, 

1.9685 

it 

11.86 

Run  15 

II 

-.2245 

lx 

.3592 

BRR 

37.4  .1132 

1.9435 

11.86 

-.2495 

.3342 

TUST 

37.4  .1257 

1.9560 

11.86 

-.2370 

.3467 

XOH 

37.4  .1374 

1.9560 

11.86 

-.2166 

.3425 

Data 

PRR 

Shifted 

for  pilot  delay  of 

18.58  .1741 

.Is 

1.6778 

10.65 

.2419 

1.2626 

BRR 

18.58  .1491 

1.6528 

10.65 

.2169 

1.2376 

TUST 

18.58  .1616 

1.6653 

10.65 

.2294 

1.2501 

XOH 

18.58  .1728 

1.6653 

10.65 

.2395 

1.2501 

Data 

PRR 

shifted 

for  pilot  dalay  of 

4.72  -1.0083 

.2s 

1.9380 

12.15 

-.1272 

.4424 

BRR 

4.72  -1.0313 

1.9130 

12.15 

-.1522 

.4214 

TUST 

4.72  -1.0188 

1.9255 

12.15 

-.1397 

.4339 

XOH 

4.72  -.9497 

1.9255 

12.15 

-.1236 

.4338 

Pilot  aodol  parameters  Identified  for  Run  2  end  Run  5  (Segment  3) 


Table  10 


Pilot  »2 
Run  91 


SlflBR&lL 

K 

II 

li 

li 

1 

-.0090 

3.3769 

-0.0267 

-0.0480 

2 

-.0316 

.1388 

-0.0266 

.1860 

3 

-.0409 

.0310 

-  .0274 

.3991 

Run  92 

Sumnt. 

It 

li 

li 

la 

1 

-.0112 

4.714 

-0.0272 

-0.0009 

2 

- 

- 

- 

- 

3 

-.0105 

0.0359 

-0.0255 

0.2718 

Pilot  Model  parameters  Identified  using  pilot 
nttdel  2  foe  Run  1  end  Run  2. 


(Tustln's  blllneer  rule) 
Table  11 


Pilot  14 
Ron  •! 


Segment  II 


K 

Ll 

T_ 

T- 

PRR 

4.0 

-.2731 

.1190 

.7503 

BRR 

4.0 

.0355 

4.1837 

-3.3644 

TU3T 

.1572 

3.3766 

.8923 

-0.0480 

ZOH 

4.0 

5.5281 

.1060 

.7377 

Segaent  92 


K 

Ll 

Iz 

Ii 

PRR 

- 

- 

- 

- 

BRR 

6.4 

-.1985 

.4126 

-.1263 

TOST 

.2573 

.1398 

.1253 

.1860 

ZOH 

- 

- 

- 

- 

Segaent  93 

K 

Il 

Il 

t3 

PRR 

9.46 

.0360 

.0923 

.4670 

BRR 

9.46 

-.4116 

.5278 

-.0185 

TOST 

.3724 

.0310 

.1352 

.3991 

ZOH 

9.46 

-.2685 

.0791 

.4544 

Pilot  Nodol  pmattri  identified  using  pilot  aodel  4  fot  Ron  1 
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Pilot  14 

Run  12 

Segment  91 

K 

li 

li 

li 

TRR 

-2.2281 

-.8781 

-.5057 

.0723 

BRR 

-2.2281 

-.0116 

4.2180 

-4.7015 

TUST 

.0877 

4.7141 

-.4576 

-0.0009 

ZOH 

-2.2281 

.6599 

.0589 

-0.5181 

Segment  93 

K 

li 

!z 

T- 

PRR 

- 

- 

- 

- 

BRR 

1.7294 

-.2843 

.0478 

-.0234 

TUST 

.0681 

.0359 

-.2223 

.2718 

ZOH 

_ 

Pilot  Model  parameters  identified  using  pilot  SK>del  4  for  Run  2 


Table  13 


Discussion  oi  Results 

The  section  will  discuss  the  identifications  of  Table 
5  through  Table  13  and  their  relationship  to  the  pilot 
ratings. 

Neal -Smith  Theory  predicts  that  as  the  aircraft  time 
delay  increases  the  pilot  must  compensate  for  this  by 
adding  more  lead.  Adding  more  lead  required  a  higher  work 
load  and  a  higher  work  load  leads  to  higher  pilot  ratings 
on  the  Cooper-Harper  rating  scale.  This  was  seen  to  be 
the  case  for  the  aircraft  flown  in  this  test.  The  higher 
the  time  delay  added  to  the  aircraft  model  the  worse  the 
pilot  ratings  became.  This  is  a  well  known  fact  in  flying 
qualities  research  and  was  expected. 

The  exact  pilot  model  parameters  predicted  by 
Neal-Smith  were  not  identified.  This  may  have  been  due 
in  part  to  the  error  bias  introduced  by  the  pointing 
accuracies  of  the  terrain  board  camera.  This  may  also 
have  been  due  to  the  fact  that  the  pilot's  attention  was 
diverted  to  the  lateral  task.  However,  if  the  pilot 
model  parameters  identified  in  run  1  and  run  2  are 
established  as  a  baseline  of  comparison  for  run  4  and  run 
5  some  Interesting  trends  can  be  seen.  Since  it  is 
not  clear  at  what  frequency  the  pilot  is  operating,  the 
lead  and  lag  will  be  examined  to  determine  trends  in  pilot 


compensation  cathec  than  calculate  the  pilot  angle  at  some 
arbitrary  frequency. 

Meal-Smith  Theory  predicts  an  increase  in  the  pilot 
compensation  to  compensate  for  aircraft  time  delay.  An 
increase  in  pilot  lead  and  a  decrease  in  pilot  lag  will 
yield  an  increase  in  the  angle  of  total  pilot 
compensation.  This  can  be  seen  to  be  the  case  when 
examining  the  pilot  model  parameters  identified  in  Table 
5  through  Table  7.  Pilot  model  1  is  used  in  each  case. 

The  aircraft  in  run  1  has  no  delay  while  the  aircraft  in 
run  4  has  a  100  milli-sec  delay.  Both  are  flown  at  the 
same  light  sequence.  If  the  values  of  lead  and  lag 
calculated  in  run  1  are  used  as  a  basis  of  comparison  for 
run  4  it  can  be  seen  that  the  pilot  lead  Increases  and  the 
pilot  lag  decreases  in  each  of  the  3  segments.  It  can 
also  be  seen  that  the  pilot  lead  increase  and  the  pilot 
lag  decreases  when  the  data  is  shifted  by  .1  seconds  and 
.2  seconds  to  make  up  for  pilot  delay.  The  fact  that  the 
pilot  will  increase  lead  compensation  and  decrease  lag 
compensation  for  an  increase  in  the  angle  of  total  pilot 
compensation  can  be  seen  in  Table  5  through  Table  7.  All 
discretization  techniques  provided  similar  results  for 
pilot  model  1.  Mo  technique  appeared  to  be  superior  to 
any  other. 

Table  8  through  Table  10  provide  a  similar  comparison 
for  run  2  and  run  5.  Wien  comparing  run  2  and  run  5  It  can 


be  seen  that  in  segment  1  and  2  the  pilot  lead  decrease 
and  the  pilot  lag  increase.  This  is  not  what  is  expected 
to  be  the  case.  However,  segment  3  yields  the  expected 
results.  This  cannot  be  explained  by  the  Neal-Smith 
Theory.  It  could  be  explained  by  an  unusually  quick  gross 
acquisition  leaving  the  pilot  an  unusually  large  amount  of 
time  to  apply  lag  compensation  for  good  steady  state 
accuracy. 

With  the  exception  of  segments  1  and  2  of  runs  2  and 
5,  Table  5  through  Table  10  show  an  Increase  in  the  total 
pilot  angle  contribution  that  is  predicted  by  Neal-Smith 
Theory. 

Table  11  presents  the  pilot  model  parameters 
identified  for  pilot  model  2.  The  discretization 
technigues  used  on  pilot  model  2  is  the  Tustin's  bilinear 
rule.  The  identifications  gave  negative  values  for  Kp 
and  tj.  For  this  reason  it  was  excluded  from  further 
consideration. 

Table  12  and  Table  13  give  the  pilot  model  parameters 
identified  for  pilot  model  3.  The  same  problem  exist  as 
for  pilot  model  2.  The  inaccuracies  in  r^,  the  pilot 
delay  exclude  the  model  from  further  consideration. 

Pilot  CfiMKAlS. 

Pilot  comments  and  Cooper-Harper  ratings  were 
recorded  during  testing.  Runs  1  and  2  were  flown  using 


the  aircraft  with  zero  time  delay.  The  pilot  describe  the 

aircraft  with  zero  time  delay  as  follows. 

"...did  not  use  roll  much,  slight  pitch  bobble 
in  fine  tracking,  it  gets  a  pilot  rating  of  3 
in  pitch  and  a  level  1  all  around..."  (12] 

The  pilot  gave  the  aircraft  with  zero  time  delay  a 

Copper-Harper  rating  of  three  which  is  from  Pigure  2, 

level  1  flying  qualities. 

After  flying  a  few  runs  with  the  aircraft  with  a  time 

delay  of  100  milli-seconds  the  pilot  was  asked  to  give  his 

ratings.  His  response  was: 

"  gross  acquisition  acceptable,  fine  tracking 
is  annoying. .. it  gets  a  pilot  rating  of  5  in 
both  axis ..."  ( 12 ) 

Prom  Pigure  (2)  it  can  be  seen  that  a  Copper-Harper  rating 

of  5  corresponds  to  level  2  flying  qualities. 

The  aircraft  with  a  time  delay  of  200  milli-seconds 

was  flown  and  received  the  following  pilot  comments. 

"...definitely  has  problems,  unpredictable, 
pilot  rating  of  7  in  pitch,  same  problem  as 
before,  small  bobble  not  really  a  pilot 
induced  oscillation. . ."(121 

A  Cooper-Harper  rating  of  7  corresponds  to  level  3  flying 
qualities. 

One  may  conclude  from  the  pilot  comments  that  the 
aircraft  flown  in  run  1  and  run  2  is  a  level  1  aircraft. 
The  aircraft  flown  in  run  4  is  a  level  2  aircraft  and  the 
aircraft  flown  in  run  5  is  a  level  3  aircraft. 
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VII.  Conclusions  and  Recommendations 

Conclusions 

Neal-Smith  Theory  accurately  predicts  the  increase  in 
pilot  condensation  for  aircraft  of  Increasing  tine  delay. 
This  was  seen  by  comparison  of  run  1  to  run  4  and  run  2  to 
run  5.  Pilot  model  1  provided  the  most  accurate  results. 
Modeling  of  the  pilot  time  delay  proved  to  be  difficult 
and  unnecessary.  Inaccuracies  in  the  pilot  delay 
identified  will  throw  off  the  other  terms  in  the 
identification.  Shifting  the  data  is  the  best  way  to 
eliminate  pilot  delay  from  the  data.  For  low  order 
transfer  functions  all  methods  of  discretization  provided 
similar  results  although  the  frequency  domain  plots 
indicate  that  Tustln's  bilinear  rule  was  the  best. 
Therefore,  the  simplest  method  could  be  used  without  any 
loss  in  accuracy. 

Rcctm cnflflUoog 

One  well  known  method  of  discretization  pole-zero 
mapping  was  not  attempted  in  this  work.  The  pole-zero 
mapping  technique  maps  a  pole  or  zero  in  the  S-plane 
directly  into  the  Z-plane.  This  method  could  be  used  and 
compared  to  the  results  obtained  by  the  other  methods. 

The  data  used  for  the  identifications  in  this  work 


was  not  obtained  for  a  purely  longitudinal  task.  The 


pilots  attention  was  Inevitably  diverted  froa  tiae  to  tiae 
on  alniaixing  lateral  pipper  errors.  The  longitudinal 
pipper  errors  are  never  coapletely  xero.  This  was  not  taken 
into  account  in  this  work.  Lateral  pilot  dynaaics  should 
be  taken  into  account.  There  is  a  need  to  deteraine  when 
the  pilots  attention  becoaes  diverted  to  the  lateral  task. 

Since  the  huaan  pilot  is  inconsistent  a  statistical 
approach  could  be  taken  were  the  paraaeters  are  averaged 
over  aany  runs  where  the  pilot  is  flying  the  saae  task. 

It  would  also  be  helpful  to  know  the  sensitivity  of 
the  RLS  algor itha  to  noise  and  pilot  reanants  when  trying 
to  identify  the  pilot  delay,  k  closed  loop  siaulation 
could  be  used  to  generate  synthetic  data  for  this  purpose. 

Matrixx  cannot  Identify  an  laproper  transfer 
function.  This  aade  the  use  of  a  lead  only  pilot  aodel 
laposslble.  Further  investigation  into  the  lead  only 
pilot  aodel  could  be  useful. 
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8IMULATIQH  PROGRAM  IfiBHX 
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PRCGRAfc  IDEM 

01MEN3I0N  THsT AC (0 : 1100),EC-1U 100 ),DELPC-2J 1100) *THETAC-3S 1100). 
,  EBIASCOJUOO).  DELPRCO:1100),PSI(3.1100) 

LOGICAL  FIRST 
INITIALIZATION 

EC-1)*  0. 

DfcLPC-2)*0. 

DELP<-1)*0. 

TncTAC-3)*C. 

TmETAC-2)*C. 

THcTA(-l)=C. 

Du  30  J  =  1.2 
DJ  3 C  1*1.1100 
PSICZtJ)®). 

CONTINUE 

DC  1C  1*0. HOC 

c  Cl) *  0. 

E&IASCI )®Q • 

DcLPCI)=U. 

DELPR  C I ) =0 . 

TnEl ACI)®C. 

TrlETACCI)®C. 

CuNTliJCr 


t  ° 

jt 


Scl  INPUT  VA4IAEL*  TIvE  PI  STORY  C  RSFEfiFNCE  NEAL-SPITH  PG  23  ) 


« c 


DJ  11  J  = 
IFCJ.GE. 
IFCJ.GE. 

lr(J«gCa 

IFCJ.GE. 
It( J«o>» 
1FCJ.U». 
IFCj.oC. 

I  r  (  j  •  u  _  • 

IFCJ.GE. 
IFCJ.Gs • 
IFCJ.'jL. 
IFCJ.Gt. 
IFCJ.GE. 
IFCJ.Gt. 
IFCJ.GE. 

lr( J.g«. 

IFCJ.GE • 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
IFCJ.GE. 
CONTINUE 


C  ,  1 1  0  u 
00  •  AND  , 
EC  .AM, 
ICC  .AND, 
1 7  C  •  AM  , 
190  .AM, 
23  C  .AM, 

3  0  C  .AND, 
34  C  .AND, 
L1C  .AND, 
*2C  .AM, 

4  4  C  .AND, 
see  .AND, 
56  C  .AND, 

5  »C  .AND , 
61 C  .AND, 
63  C  • ANC , 
72C  .AND. 
77C  .AND, 
880  .AND, 
91 C  .AND, 
960  .AND, 
980  .AND, 
10C0.ANC, 
10X0. AND, 
1020) 


J.LT.30 

) 

THETACC J)* 

0.0 

J.LT.10C 

) 

TPETACCJ)* 

-1.0 

J.LT.1T0 

) 

THETACCJ)® 

-2.0 

J.LT.19C 

) 

THETACCJ)* 

0.0 

J.LT. 239 

) 

THETACCJ)* 

-1.0 

J .LT .300 

) 

THETACCJ)® 

0.0 

J.LT. 340 

) 

THETACCJ)® 

4.0 

j.lt.lio 

) 

THETACCJ)® 

1.5 

J.LT.A20 

) 

THETACC J)* 

-1.0 

J.LT. 440 

) 

THETACCJ)® 

-1.5 

J.LT.500 

) 

THETACC J)* 

-2.0 

J.LT. 560 

) 

THETACCJ)® 

0.0 

J.LT . e  90 

) 

THFTACC J)» 

-2.5 

J.LT. 410 

) 

THFTACC J)» 

0.0 

J.LT. 630 

) 

THETACCJ)® 

3.0 

J.LT. 720 

) 

THETACC J)* 

0.0 

J.LT. 770 

) 

THETACCJ)® 

5.0 

J.LT. 830 

) 

THETACC J)« 

3.0 

J.LT. 910 

) 

THETACC  J)* 

4.0 

J.LT. 960 

) 

THETACCJ)® 

-1.0 

J.LT. 930 

) 

THETAC  C J)* 

-1.5 

J.LT. 1000) 

THET  ACC J)® 

-2.0 

J.LT. 1010) 

TMETACCJ)* 

-1.0 

J.LT. 1020) 

THETACC J)* 
THETACCJ)* 

-0.5 

0.0 

CALCULATE  CCEFFICIENTS  FOR  THE  OPTIMUM  PILOT  MODEL 
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c 

c 

<i> 

c 

c 

< 


c 


c 


TAt*  .1 
TALLAG  s  .3 
TALI s* 1933 
TAt2*.01 
T  AL3  *T  AULAG/2 • 

RKP  =  7 • 47  26 
RKP*. 2942 

OENOH  *  <  T  A  U2  *  T  AU  3  ) 

AG  *  RKP/OENDM 

A1  *  RKP  *  <  TAbl  -  TAt'3  )  /  DFNCM 
A 2  =  -RKP  *  C  TAU1  *  TAU’  )  /  r§NCM 
o  C  =  l./CsKwf* 

bl  :  (  T Ati2  ♦  TA'J?  )  /  CENCM 

03K3V2  =  C  1./<TAU*TA'J)  ♦  Fl/TAU  *  90  > 

A  AO  -  (  A2/(TAU*TA'J)  ♦  Al/TA  U  ♦  AO  )/0EN3M2 
A A1 =  C  -2.*aZ/(7Au*TAJ)  -  AJ/TAU  3/CENCM2 
A  A2=  (  A2/(T  AUvTAU)  3/DEN3»»2 
ofc  1*  (  2./CTAL'#TAU>  ♦  nl/TAL  J/DFNCM2 
862=  C-l  ./(TAU*TAo)  )/:e\'0*2 


PRIM 

CUT  CC= 

PP2CIEWTS 

F  3C 

PILOT  MODEL 

PH;,T 

v ,  •  • 

PRIM 

*»  '  T  £  L 

=  •  ,  T  £  *J 

prim 

*.  'TAU1 

=  •  »  T  £  U 1 

PRIM 

*  »  '  T  A  o  2 

=  '  »  T  AU2 

PRIM 

*  *  *  T  £  'J  3 

=  * i T  AU3 

Phi '« T 

v,  'RKP 

*  •  ,  RK? 

PRIM 

*,  'AC 

*  *  »  A  C 

PRINT 

*,  *  A1 

=  '  »  £  1 

PRIM 

v t  '  A2 

=  1  »  A  2 

PRIM 

*  «  '  i  C 

* 

PRINT 

*»  '£  1 

*  '  »E1 

PRIM 

V|  *  M  *•  u 

=  •  ,  A  A  0 

PRINT 

'  A£1 

*  • , AA1 

PRINT 

'  A  A  2 

=  •  »  A  A  2 

PRINT 

$i  *361 

*  *  f  ?M 

PRINT 

2 

=  *  till 

PRINT 

*,  i  t 

CALCULATE  CGE 

FcICIcNTS 

FJ5 

TM?  AIRCRAFT  M3CEL 

RM.*T1nCH*.C3937 

RlNCnTl*'^.*2i*<* 


CG*0.0 
Cl*  9*0 
C2*4.2 

00*1  *1376**  INC  MT  MM. 
Q1*.948**IKCHT*M 


m 

« 


DENCH*  Cl./TAo**3.  ♦  C2/TAL##2.  ♦  Cl/TAU  ♦  CO  > 
CC1*  -<  -3./TAL*$3.  -2.*C2/TAU**2.  -Cl/TAU  >/DENOM 
CC2*  -<  3./T  AU**>3.  ♦  C2/TAU**2.  7/OENCP 
CCS*  <  1./TAUPP3.  >/D«ND* 


0 


129- 


000*  C  OI/TAU  ♦  DO  J/OSNCM 
001*  C  -Jl/T-0  >/*ENCP 


PKInT 

CUT  CC  = 

FFICItNTS  PC?  TH£ 

AIRCRAFT 

PRINT 

«•  *  1 

PRINT 

*,*cc 

S 

SCO 

PRINT 

*.  *C1 

* 

SCI 

PRINT 

#t'C2 

9 

SC2 

PRINT 

*.  *oc 

X 

•  t  ro 

PRINT 

♦  *  •  01 

9 

sci 

PnlNT 

*t'CCl 

X 

scci 

PRINT 

*,  *CC2 

X 

SCC2 

print 

*•  'CCA 

X 

sec? 

PRINT 

*»  »ocj 

S 

SCI-0 

PRINT 

X 

SCL1 

PRINT 

9,  •  • 

QPtN 

file:  -c 

*  ^ 

ATA  STCPAC-E 

OPcNC 

2,  PILE** 

.STATUS*' NEW) 

OPENC 

3  »  r IL  E  *  ' 

2?L 

r . • , STATUS* • NEW 

'  ) 

. . .  'COTCAT.' ,STATUS*'N?W) 

WRIT  c(4»lul2) 

WRITt(4,1013) 

WkiT:(<.»101fc) 

1012  PuRKATClX,'  ') 

1  Cl  3  FuR^TliX,'  S?X«’C?LC(IN)  S5X,*  THETA  (DEG) 

:  OPERATE  LCC? 

rikiT  s.TPL: • 

DO  2  0  K « 1 ( 1 1 0 0 

k  cCKJ=Th;1 ACCiO-Tt-FTACOl) 

CcLP(N)*;::* Ci LP(*-2>+?: l*:flP(K-l) 

♦J-2vr(k-2)^Aaiwr(K-l)46A0«ECK) 

0*D:L=C<-1 ) 

C-LL  ‘-liCj.t 
THETA(K)*Y 

F  *  1  N  T  *  r  ' 

WRIT; 

ruRHATOS 
HXTsC2»l0lG3r<O 
W*lTE(j» lC10)ufLc(N) 

PO*:',AT(;  lfc.O 
CONTINUE 


» 5X) 


»  P I  *  ST  ) 


(*)*Y 

*,'  =  *  f^.LP*  »?HET  A*  *  ,E(X),CELPClOtTHETA(IO,THETAC(K) 

CotlOll):  OOt?ELP(*)tTH*T«CK>vTHETACCK> 

.TClX*61^.4(SX*c12.*«5Xt«12.4v5XfP12.S> 

C2(lUl03rCO 


ma&Mmam&imM 


ATC2»1)*0. 
ATC2,2)»-.$6l3 
AT( 2 , 3  )*  •  0 £  Q-* 
AT(3,1)=C. 

ATC3.25  —  .7240 
AT(3  »  53* • 62  3  5 

5T(1)«. 00016777 
dT(2)».0042000G 
6T(3)*«08040C0C 

CTCl )«1. 1376  *25.* 
CT(2>«.94bC  *  23  •* 
CTC33=  C. 

XTu(l)=j. 

XTCC20-3. 

XTCC3)=0. 

XTC13=G. 

XT(2)=C. 

XTti 3=0. 

E?4i  ir 

0  0  ir  1  =  1,3 
X7(l)s0. 

25  CLMIU'w- 

U0  1C  1=1,3 
DC  15  J*  1 , 3 

T  ;  I  -  .  r  =  -  7  C  1 ,  J  )  -•  •  X  T  C  C  J  ) 
XT  C  I )  *X7  C  I  )  +  T  r f  ? 

C  CM  1 1,0  r 

XT(:)  =  x7C)4r  T(l)»  j 
CQM  I/mLc 

Y  *0  . 

DO  20  1=1,3 

xtuci)=xtc:) 

YT  =  CT ( 1)$X7 ( I  ) 

Y  =  Y  4  Y  T 

20  COM  If,  Ur 

RETURN 
C  .Mi 


<9 


cRkCkCOEu  ) 


l=lpc:n> 


THETft(DEG) 


THETA c  (p£<*) 


0.0330 

0.0000 

0.0000 

C.0000 

0.0000 

O.COOO 

o.cooo 

C.00C3 

0.0300 

0.0000 

0.0000 

0.03 0  0 

0  .  vj  0  0  0 

C  .  C  0  0  0 

0.0000 

o.ocoo 

0.3000 

C  .  0  0  0  c 

0.0030 

0.0000 

0 . 0  J  0  0 

0.0003 

0 . 0 C  vO 

C .0  ooo 

C.OuOC 

0.0000 

0 . 5  3  C  0 

o.cooo 

o .  c  o  o : 

C.0032 

o .  c :  o : 

0.0000 

o.cooo 

0.0300 

0  •  0  0  J  0 

o.cooo 

3  .  C  u  3  w 

0 .0003 

0  .  C  3  0  0 

0.0000 

0  .  C  o  u  v 

C  .  0  3  0  0 

0.0030 

0 . 0  0  0  c 

0.0330 

o.cooo 

o.occo 

C.0000 

0 . 3  3  w  3 

•3.0  30  0 

o.ocoo 

0.0000 

W  .  u  w  o  w 

o.cooo 
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0.2075 

3.3189 

4.00000 

0.6811 

0.2C40 

3.3301 

4.00000 

0.6699 

0.2909 

3.3406 

4.00000 

0.6594 

0.1975 

3.3514  * 

4.00000 

0.6486 

0.1945 

3.3615 

4.00000 

4.5615 

0.9757 

3.3719 

-1.00000 

<*.3719 

-1.3536 

3.4663 

-1.00000 

4  .4  o  6  o 

-1.445* 

3.2897 

-1.00000 

4.2357 

-1 .4601 

3.1562 

-1.00000 

4.1562 

-1.354? 

3.0301 

-1.00000 

4.C3  Ji 

-1.2795 

2.9430 

-1.00000 

3.5430 

-1.215° 

2.8572 

-1.00000 

i.€=7fc 

-1.17:3 

2.7921 

-1.00000 

5.7*21 

-1.1442 

2.7199 

-1.00000 

3.71** 

-1.1230 

2.66C7 

-1.00000 

3  .  t  o  C  7 

-1.0964 

2.5933 

-1.00000 

3.53 3i 

-1.0515 

2.5361 

-1.00000 

3.53ol 

-1.039? 

2.4716 

-1.00000 

5.4716 

-1.043? 

2.4158 

-1.00000 

3.4153 

-1.0236 

2.3537 

-1.00000 

3.3i27 

-1 .00*1 

2.2994 

-1.00000 

3.2394 

-0.9?5? 

2.2398 

-1.00000 

3.43*2 

-0. 97.-4 

2.1869 

-1.00000 

3.1665 

-0.5552 

2.1297 

-1.00000 

3.12*7 

-0.5402 

2.0783 

-1.00000 

3.0732 

-0.9227 

2.0232 

-1.00000 

3.0232 

-0.5  0t2 

1.9733 

-1.00000 

2.5733 

-0.8914 

1.9204 

-1.00000 

2.52  0** 

-2.6772 

1.8720 

-1.00000 

2.5720 

-  0 . 6  c  1 1 

1.3211 

-1.00000 

2.6411 

-  0  •  fc  4  7  3 

1.7742 

-1.00000 

2.7742 

-0 .3315 

1.7252 

-1.00000 

2.7224 

—  0  .  S  1  :  4 

1.67S7 

-1.00000 

2.6797 

-C.8036 

1.6325 

-1.00000 

2 .6323 

—0.7505 

1.5S84 

-1.00000 

2 • 5  0  34 

-0.776? 

1.5430 

-1.00000 

2.5430 

-  0 . 7  o  3  f 

1.59C3 

-1.00000 

2.5003 

-C .7404 

1.4566 

-1.00000 

2.430b 

-0.7376 

1.4152 

-1.00000 

2.4152 

-0.7244 

1.3731 

-1.00000 

2.37 jl 

-0.7124 

1.3330 

-1.00000 

2.3330 

-0.699? 

1.2924 

-1.00000 

2.292* 

-0.68  32 

1.2537 

-1.00000 

2.2537 

-0.6760 

1.2145 

-1.00000 

2.2145 

-0.6648 

1.1770 

-1.00000 

2.1770 

-0.6531 

1.1392 

-1.00000 

2.1392 

-0.6421 

1.1030 

-1.00000 

2.1030 

-0.6309 

1.0665 

-1.00000 

2.0665 

-0.6203 

1.0315 

-1.00000 

2.0315 

-0.6094 

0.9963 

-1.00000 

1.9963 

-0.5992 

0.9524 

-1.00000 

1.3624 

-0.5587 

0.9285 

-1.00000 

1.9265 

-0.573? 

0.8958 

-1.00000 

-148- 


-1.8958 

-0.56S7 

0.8630 

-1.00000 

-1.8630 

-C.5592 

0.8313 

-1.00000 

-2.3313 

-C.4710 

0.7998 

-1.50000 

-2.2996 

-0.6945 

0.7776 

-1.50000 

-2.2776 

-0.6939 

0.7285 

-1.50000 

-2.2263 

-0.6861 

0.6846 

-1.50000 

-2.1346 

-0.6662 

0.6416 

-1.50000 

-2.1416 

-0.6497 

0.6033 

-1.50000 

-2.1033 

-0.6343 

0. 5654 

-1.50000 

-2.0654 

-  U  .  6  2 1  9 

0.5304 

-1.50000 

-2.0304 

-C .6?98 

0.4949 

-1.50000 

-1.S949 

-0.5393 

0.4614 

-1.50000 

-1.9614 

-0.5 3 “5 

0.4273 

-1.50000 

-1.9273 

-C.5757 

0.3949 

-1.50000 

-1.6349 

-0.5694 

0.3621 

-1.50000 

-1 .6c21 

-G . : 5  9  0 

0.3303 

-1.50000 

-1.6306 

-0.5491 

0.2991 

-1.50000 

-1.7931 

-0.5400 

C.2683 

-1.50000 

-1 • 76  o e 

-0 .6306 

0.2382 

-1.50000 

-1.72:2 

-C  .5217 

0. 2089 

-1.50000 

-1.73*9 

-0.5126 

0.1795 

-1.50000 

-1.679: 

-0 .50^1 

C.  1511 

-1.50000 

-2.1311 

-3.4145 

0.1227 

-2.00000 

-2.12^7 

-0.6414 

0.1038 

-2.00000 

-2.103= 

-0 .6416 

0.0576 

-2.00000 

-  c  •  0  a  7 1 

-0.62-43 

o.oiee 

-2.00000 

-  2  .  C 1 6  3 

-C.615S 

-0.0232 

-2.00000 

-1.976: 

-  C  .  -5  0  C  2 

-0.0585 

-2.00000 

-1.9413 

-  c  .  5  c  5  7 

-0.0935 

-2.00000 

- 1  •  9  0  o  2 

-0.574: 

-0.1257 

-2.00000 

-1.67-3 

-0.5623 

-0.1585 

-2.00000 

-1  .  c -  1  3 

-0.633? 

-0.1392 

-2.00000 

-1.8103 

-0.5432 

-0.22C6 

-2.00000 

-1.77*4 

-  0  •  5  5  4  3 

-0.25C3 

-2.00000 

-1.7497 

-0.5243 

-0.2305 

-2.00000 

-1.719; 

-0.31:2 

-0.3093 

-2.00000 

-1.6  907 

-0.5071 

-0.3385 

-2.00000 

-1.6615 

-0.4337 

-0.3663 

-2.00000 

-1.6237 

-0.4  300 

-0.3944 

-2.00000 

-1.60:6 

-0.4-19 

-0.4213 

-2.00000 

-1.5737 

-0.4775 

-0.4485 

-2.00000 

-1.5315 

-0.-556 

-0.4745 

-2.00000 

-0.5255 

-0.6144 

-0.5007 

-1.00000 

-0.4*93 

-0.1412 

-0.5429 

-1.00000 

-0.4571 

-C  .11 59 

-0.53C7 

-1.00000 

—  0  •  4  o  9  2 

-0.106? 

-0.5264 

-1.00000 

-0.4736 

-0.1204 

-0.5237 

-1.00000 

-0.4763 

-0.1290 

-0.5279 

-1.00000 

-0.4721 

-0.1350 

-0.5324 

-1.00000 

-0.467o 

-0.1361 

-0.5404 

-1.00000 

-0.4596 

-C .1366 

-0.5468 

-1.00000 

-0.4532 

-0.1347 

-0.5553 

-1.00000 

0.0553 

-0.211? 

-C.5619 

-0.50000 

0.0619 

0.0236 

-0.5786 

-C. 50000 

0.0786 

0.0339 

-0.5664 

-0.50000 

0.0664 

0.0376 

-0.5599 

-C. 50000 

0.0599 

0.0265 

-0.5527 

-0.50000 

0.0527 

0.0232 

-0.5505 

-0.50000 

0.0505 

C  .0102 

-0.5472 

-0.50000 

0.0472 

0.0165 

-0.5469 

-0.50000 

0.0469 

0  •  C 1  4  5 

-0.5448 

-0.50000 

0  •  0<*4d 

0.0144 

-0.5448 

-0.50000 

0.6448 

-0.0652 

-0.5431 

0.00000 

0.5431 

0.1675 

-0.5516 

0.00000 

0.5516 

0.1759 

-0.5313 

0.00000 

0.5315 

0.1771 

-0.5169 

0.00000 

0.5169 

0.1656 

-0.5019 

0.00000 

0.5019 

C.1579 

-0.4921 

0.00000 

0.4921 

0.1507 

-0.4812 

0.00000 

0.4812 

0.1469 

-0.4736 

C. 00000 

0.4736 

0.1425 

-0.4642 

0.00000 

0.4642 

0.1402 

-0.4571 

0.00000 

0.4571 

0.1370 

-0.4484 

0.00000 

0.4424 

0.1351 

-0.4415 

C. 00000 

0.4415 

C.1322 

-0.4332 

C. 00000 

0.43 32 

0.1304 

-0.4264 

0.00000 

0.4264 

0.1277 

-0.4194 

0.00000 

0.4134 

0.1259 

-0.4118 

C. 00000 

0.4113 

0.1233 

-0.4042 

C. 00000 

0.404^ 

u  •  1  4  i  > 

-0.3977 

c.oooco 

0.3277 

0.1171 

-0.3904 

0.00000 

0.3904 

0.117*. 

-0.3841 

0.00000 

0.3341 

0.1151 

-0.3771 

0.00000 

0.2771 

0.1132 

-0.3709 

C. 00000 

0 . 3  7  C  9 

0.1112 

-0.3642 

0.00000 

0.3642 

0.1095 

-0.3582 

0.00000 

0.3562 

0.1074 

-0.3518 

0.00000 

0.351c 

0.1C67 

-0.3460 

C. 00000 

0.3460 

0.1037 

-0.3396 

0.00000 

0.335c 

0.1021 

-0.3341 

0.00000 

0.3341 

0.1002 

-0.3292 

C. 00000 

0.22  52 

0.0996 

-0.3227 

C. 00000 

0.3227 

0  .  C  9  6  9 

-0.3170 

0.00000 

0.3170 

0.0552 

-0.3117 

C. 00000 

0.3117 

C  .  0  5  3  5 

-0.3061 

C. 00000 

0  •  3  0  c  1 

C  •  J  *  7 

-0.3010 

0.00000 

0.301Q 

C  .  0  5  0  3 

-0.2957 

0.00000 

0.2*57 

0.093? 

-0.2907 

0.00000 

0.2307 

0 . 0  w  7  2 

-0.2956 

0.00000 

0 . 2  o  2  c 

0 . 0  5  5  r 

-0.28C8 

C. 00000 

0.2308 

0.0942 

-0.275S 

0.00000 

0.2755 

0.0525 

-0.2712 

0.00000 

0.2712 

0.0513 

-0.2664 

0.00000 

C.2064 

0.0500 

-0.2619 

C. 00000 

0.2619 

0.07*6 

-0.2573 

0.00000 

0.25  73 

0.0773 

-0.2530 

0.00000 

0.2530 

0.0759 

-0.2495 

C. 00000 

0.246a 

0.0746 

-0.2443 

0.00000 

0.2443 

0.073? 

-0.2400 

0.00000 

0.2400 

0.0721 

-0.2360 

0.00000 

0.2360 

0.0709 

-0.231? 

C. 00000 

0.2318 

0.0696 

-0.2279 

C. 00000 

0.2279 

0.0684 

-0.2239 

0.00000 

0.2239 

0.0672 

-0.22C1 

0.00000 

0.2201 

0.0660 

-0.2163 

0.00000 

0.2163 

0 .0649 

-0.2126 

0.00000 

0.2126 

0.0638 

-0.2089 

0.00000 

0.2089 

0.0627 

-0.2053 

0.00000 

0.2053 

0.0616 

-0.2017 

0.00000 

0.2017 

0.0606 

-0.1983 

0.00000 

( 

0.1383 

0.0535 

-0.1948 

0.00000 

0.1348 

0.0585 

-0.1915 

C. 00000 

c 

0.1915 

C  .  0  5  7  5 

-0.1882 

0.00000 

0 • 1 8  £2 

0.0555 

-0.1350 

0.00000 

0.1350 

0.0555 

-0.1818 

0.00000 

(§ 

0.181s 

0.0546 

-0.1766 

0.00000 

0. 1 7  ofc 

0.053* 

-0.1755 

C. 00000 

0.1755 

0.0  =  27 

-0.1725 

C. 00000 

« 

0.1725 

0.0=18 

-0.1696 

0.00000 

0.1696 

0 . 0  5  0  3 

-0.1666 

0.00000 

0.1666 

0.0500 

-0.1638 

0.00000 

c 

0.1t33 

Ij  •  J  **  7  1 

-0.1609 

C. 00000 

0.1609 

0.04V? 

-0.1592 

0.00000 

0  •  1  v  c  <» 

0.0475 

-0.1554 

C. 00000 

c 

0.1554 

0.0465 

-0.1528 

C. 00000 

0.1525 

C  .  0  4  5  8 

-C.1501 

C. 00000 

0.1501 

0  .  C  4  5  1 

-0.1475 

C. 00000 

< 

0.1475 

0.1-50 

0 . 0  4  4  3 
0.0^35 

-0.1450 

-C.1425 

0.00000 

0.00000 

O.l^Is 

0 . 0  4  2  3 

-0.14C0 

C. 00000 

r 

0  •  1  4  J  W 

C  .  0  4  2  0 

-0.1376 

C. 00000 

V 

0.157s 

0.041? 

-0.1353 

0.00000 

0.1352 

0.04  o £ 

-C.  1329 

C. 00000 

( 

< 

« 

fS&tk. 

* 

c 

c 

« 

c 

a 

a 

c 

-152- 

c 

■:  \v.  v>y*y»y*y» 

,  i  *  i '  •  1  ^  *  i  * '  * .  *  » . 1  ;• 

UV  UmiH  BASK 

IS1 


PROGRAM  BASE 


LOGICAL  FIRST, CELAY 
CHARACTERS^  ANSWER 

fWRlTSC*,201) 

FuRMATClX, 'PILOT  Nlf'EER  C  £  ♦  1  ♦  2  ♦  3 . 4  ,  ft  )  •  > 

READ  (*»2G2),aNSr£R 
202  FORHATCAd) 

4*1.4445 
8*7.4725 
C*.Qi 
J  =  1.0 
es*i 
r  =  .  1  3 


AC*  a*C 
AO*  A*3 

AC*  A V  Z 

A  p*  ASr 
d  C*  c  vC 
50*  5*0 
cc*  5*5 
3r  =  t*F 
Cc*  C*E 
OE*  0v£ 
up*  Dvr 
FC*  fv C 

p - =  rvr 
KTC£'j*37.3 

FIRST*.TRUe. 

0 EL AY*,  p ALEE. 

upen C2,. =  :le*' case. cat  a*, status* 'cl:'> 
RRITc(2,1C1) 

WRIT  5(2 , 103 } 

00  IOC  1*1,100 


301 


xo*i/io. 

?  H I  *  E  *  X  3 
rfilJ*PhI*HTucv» 

C?  *CCS( CF I  ) 

SP*5INCPFI) 

CcNOM*  C*C*XC*XC-»0*0 

IF(ANS*?S. =2. • 3' )  THEN 
1F(PIrST)*h:TE(*,301> 
FORMAT  CiX  » •iASSLINt') 


XR«CCAC-:C)*XO*SCt-<  AC*XO*XO*SC)*CP)/CENuH 
XI«C C AC- EC )* Xj *CP-C AC* X0*XC+E0 )*«?)/: ENOP 


XMGT*SCRT( Xk**2.*XI*m2.) 
XPHT*ATAN2(XI»XR) 


302 


ELSE  IFCANSWEP.EC.  'I')  TpEN 
IF<FIft$T3wRITE($,302) 
P0RMATC1X, 'PILOT  1') 


X«N*B 

XIN*A«XO 


XRD»Q 
XXD«OXO 

ELSE  I=CfiNSwcQ.£C.*2«)  THEN 
ZFCFIftST)*iRIT£(#»303) 

303  FQRMATC1X, 'PILOT  2') 

XRN*(  AE^Xu^XOSO 
XIN=(  A-3*E)»XG 

XR  0  =  0 
X±3  =  C**sXtj 

ci.SE  IFCANShieR.SC.'3*>  TPEN 
Ir(rIRST)WRITE<*,3CO 

304  rCNMATClX,  •  PILOT  3') 

X  ft  N  *  a 
X  _  N  =  A  v  X  u 


306 


A  &  “  u  ~  i-  -  v  X  j  •>*  X  0 
X1G=CC +Gc)«XC 

ELSE  IFCANiwcft.s:.'*')  T  *■  r  N 
IP<r;ftST)nWlTcC#*3C5> 

FORtfATC IX, 'PILOT  <• '  ) 

X  R  \  s  A  r  5,«  X  L  w  X  3  +  a 
XIN*( A-:f )«X j 

XR3* 3-FC#X jvXO 

X I  C  =  (C+CP  )  *XXC 

ELSE  IF(AT«6W-‘;#Si«*>J')  TPfN 
IPCFIRST)  t-en 
URITE(*,30e>3 

P0RMATC1X, 'filCHCLS  PLOT  C  AT  A ' ) 
CELAY*.TR'j3. 

Enc  Ie 


C  3PTIMUM  PILOT  VCOfL  AT  Cm£GA  =  2 


C  XiiN  =  -133. Y  «  *2*XC  ♦  3  5  0 . 1 

C  XIN  1  c  7  2 . 7  -i  Xu 


C  XftO  *  XjftX.5 XX'J#X0  -  <.23.vXC*XO 

C  XZO  ■  -10*.2*X.-*XDm,T0  ♦  2t3.?*xC 


BASELINE  AIPCPArT  A N C  TI*:  L  ELAY 

X RN  *  1.1376 
X IN  *  .'J43SXX0 


C  XRD  »  -4.2*XC*XC 

C  XIO  =  9.SX0  -  XCvXQ^XQ 

C  OPTIMUM  PILOT  MCOSL  AT  C*E&A  £  - 

fRNC  *  3:0.s 

RN1  *  iC3.2  <  , 

RN2  *  2 to.-* 

RD4  »1. 

R03  »  1C4.2 
RD2  *  429.0 
ROl  *  8S9.& 

XRN  *  kNO  -  kN2*XO*XQ 
XIN  =  kM  *XG 

XRD  "  XC~51S4  ~  C 2 *■' X C ^ 2 
XIO  =  RC1*XC  -  R02«XG#*2 

t NO  IF 

XMGTN  =  i-.RTC  XCAvis2.  +  x:n*«2) 
X?*lTN  =  aTi"i2<X:N,Xk\) 

X  ?  HNC  =  XP.-TN*RTCEG 

XMGTC  =  SjRTC  X*C*«2.-fXI3S*2) 
X^-TC  =  ATu,?(XlC»X?C) 

XPhJG  =  X  r 4 T  0  » R  T  0  E  G 


XMGT*  X:'G7  -i/X.-ITC 
X-HT*  X  ®  F  T  N  -  X  P  PTO 
IrCC.LflY)  X5hT»XCMT-PHl 


401  CONTINUE 

IFCCISST3  =:^ST=.FSL$r. 

X  /  G  =  2  0  .  #  L  'J  G  1  0  C  X ."  C  T  ) 

XPr*X=hT*KTOfG 
IFC-lNi^r  5.  .‘N*  )  T  ^  c  N 

IFCXPM.jT.-T.)  XPp=XPF-2cO. 

END  I  r 

«k  IT  E c  2  » 1 0  5  )  X  0  *  X y G *  X pH,  X s  N  t  X I N  *  X  P^NC  ♦  X »C  *  XI  J ,  XP**0 D  1  P h  1 0 
W  R  I T  c  £  *  » 1 0  2  )  X  0  *  X  * -1 G  ♦  X  P  M 

100  CONTINUE 
CL J  Sc  C  2 ) 

101  rCiRMAT<7X,*  C-'cGA  •  ,5X ,  •  MftGCOi) '  «2X«  '  pPASE  C0:G)  %2Xf»  XRN  », 

.  2  X ,  •  X I N  NLP  ’  ,  2  X  t '  ‘  =  0  *fiX,'  X I D  'f2X,'  -UN  %2X,'  PHI  ') 

FORMAT  C  IX  »  i  F  1  .4) 

FCRMATClXf  10P10.2) 

FORMAT (IX ) 


102 

105 

103 


fMQUMICV  BEflPOMflE  OF  OPTIMAL  PILOT  MODEL 
UfllMC  THE  ACTUAL  TIME  DELAY 


OMEGA 


MAGNITUDE  (OB)  PHASE  (CEG) 


0.0000 

0.1000 

0.2000 

0.3000 

0.4000 

0.9000 

O.BOOO 

0.7000 

O.BOOO 

0.9000 

1.0000 

1.1000 

1.2000 

1.3000' 

1.4000 

1.5000 

1.6000 

1.7000 

1.8000 

1.9000 
2.0000 
2.1000 
2.2000 

2.3000 

2.4000 

2.5000 

2.6000 

2.7000 

2.8000 

2.9000 
3.0000 

3.1000 

3.2000 

3.3000 

3.4000 

3.5000 

3.6000 

3.7000 

3.8000 

3.9000 
4.0000 

4.1000 

4.2000 
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